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A classical field theory is proposed for the electric current and the electromagnetic field interpo- 
lating between microscopic and macroscopic domains. It represents a generalization of the density 
functional for the dynamics of the current and the electromagnetic field in the quantum side of 
the crossover and reproduces standard classical electrodynamics on the other side. The effective 
action derived in the closed time path formalism and the equations of motion follow from the vari- 
ational principle. The polarization of the Dirac-see can be taken into account in the quadratic 
approximation of the action by the introduction of the deplacement field strengths as in conven- 
tional classical electrodynamics. Decoherence appears naturally as a simple one-loop effect in this 
formalism. It is argued that the radiation time arrow is generated from the quantum boundary con- 
ditions in time by decoherence at the quantum-classical crossover and the Abraham-Lorentz force 
arises from the accelerating charge or from other charges in the macroscopic or the microscopic 
side, respectively. The functional form of quantum renormalization group, the generalization of the 
renormalization group method for the density matrix, is proposed to follow the scale dependence 
bX}' through the quantum-classical crossover in a systematical manner. 
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I. INTRODUCTION 



Classical systems traditionally serve as starting points for quantization. But the opposite order, namely the con- 
struction of macroscopic physics from quantum principle, is needed to understand the great division line in Physics, 
' the quantum-classical crossover. One expects no surprise from the traditional way of thinking because the usual 
, equations of motions are supposed to be recovered in the narrow wave-packet limits according to Ehrenfest theorem. 
Thought formally correct, this theorem does not guarantee that the original, naive scenario of classical particles or 
bodies is recovered because the relativistic and many-body aspects of quantum physics introduce correlations in the 
' narrow wave-packet limit which render the picture of classically localized particles and bodies more involved. The 
I goal of this paper is to derive classical dynamics, namely action and variational equations of motion for excitations 
Qh. in QED. 

^ ' The dynamics of expectation values can be obtained by means of the closed time path (CTP) formalism proposed 
by Schwinger long time ago [1]. One expects that the dynamics for the electric current and the electromagnetic field 
will be governed by the action 



^ ; S = -^A'h I ds-^e, I A^{x)dx^ - - I d^'xidM^) - duA^ix))' + Sq (1) 



involving the world lines Xi{s) of charge e^, mass Ali and the vector potential Afj_{x) and Sq standing for the corrections 
due to quantum fluctuations. But there is a negative result conjecturing that no variational principle is available for the 
derivation of the equations of motion [2] and in fact, up to our knowledge no positive result had been communicated 
in this direction ever since. It will be shown below that the action (1) together with its variational principle can 
actually be derived. The term Sq in Eq. (1) represents the systematical improvement of classical electrodynamics 
by taking into account quantum effects arising from the polarization of the Dirac-see without interaction, vacuum 
polarizations due to electromagnetic interactions and boundary conditions in time. 

Once the action governing the expectation values is constructed a number of interesting questions open up. The 
first set of questions concerns the form of the action. How can we separate off the many-body aspects from a particle 
propagating in the non-interacting Dirac-see and what corrections will be added to the free mechanical action, the 
fist term on the right hand side of Eq. (1)? The dynamics of the expectation values of local operators should give a 
better insight into the quantum-classical transition regime because the space-time resolution of the expectation value 
of a local operator is limited by the UV cutoff only. How can we separate microscopic and macroscopic effects in the 
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expectation values, or preferably in the action? The effects of the polarization of the Dirac-see are included in the 
quantum corrections to the action. Does this happen in the usual manner polarization is taken into account in classical 
electrodynamics? Another set of questions arise about the radiational time arrow. The boundary conditions in time 
are realized on the microscopic level in QED and appear as infinitesimal imaginary contributions in Sq. How does 
the experimentally well established retarded Green function arise from the time reversal invariant dynamics? Finally, 
such a calculation addressc;s the measurement theory of Quantum Mechanics. In fact, in analyzing the measurement 
process one always ends up with the study of the interactions between a small and a large system. Are there some 
peculiarities in the dynamics when the large parameter, the ratio of the number of degrees of freedom of the large and 
the small system, diverges, like spontaneous symmetry breaking occuring in some models in the thermodynamical 
limit? We attempt to give below at least some indications about the answers. 

The CTP formalism has originally been introduced in Quantum Field Theory [1-12] for the description of the 
time dependence of expectation values. Similar scheme was developed independently in the framework of Quantum 
Mechanics, too. The formal invariance of transition probabilities under time inversion lead to the construction of time 
symmetrical Quantum Mechanics [13, 14]. The explicit appearance of the density matrix in the formalism explains 
the applications of this scheme in calculating the reduced density matrix in Quantum Mechanics [15] with special 
attention paid to dissipative processes [16-20]. Furthermore, the description of decoherence [21-24] in the framework 
of consistent histories [25-27] can be achieved in a natural manner [28-30] in the CTP formalism. Finally, one should 
mention promising applications of this scheme to cosmological problems [31-34], to kinetic Quantum Field Theory 
[35-38] and to the renormalization group [39, 40]. 

In these applications of the CTP formalism the degrees of freedom are doubled, due to the systematic implementation 
of time evolution in the Heisenberg picture and/or to the canonical description of dissipative processes. But this 
doubling remains rather formal from the point of view of classical physics. The present work is based on the observation 
that one can actually construct variational principle and canonical dynamics for the expectation values of observables 
within the CTP formalism. The choice of the microscopic degrees of freedom which become classical in the macroscopic 
region is unique up to a parameter which is allowed due to the unitarity of the time evolution for closed systems. 
The other, independent combination of the doubled degrees of freedom takes care of the coupling of the system 
to its environment and becomes suppressed in the macroscopic region. Such an explicit splitting of the variables 
offers a new, more natural starting point for the exploration of the quantum-classical crossover, presented in this 
paper within the framework of QED. The scope of this work is restricted to general, qualitative issues, such as the 
current dynamics in the Dirac-see, the polarizability of the vacuum, the radiation time arrow, the decoherence of non- 
relativistic charges and the renormalization group scheme for the density matrix to describe the quantum-classical 
crossover. The classical dynamics derived in these contexts yields new results already at such a qualitative level. We 
plan to present the detailed, quantitative studies of these and related issues in subsequent publications. 

The organisation of this paper is the following. The action for the expectation values and some partial answer to 
these questions will be found within the CTP formalism, introduced in Section II with special attention paid to the 
time evolution of the averages of observables and the role of the boundary conditions in time. The simplest system 
where the role of the boundary conditions in time can be traced in forming the time arrow is that of non-interacting 
particles. Section III contains the derivation of the action for the expectation value of the free photon field. The 
electric current is a composite operator and its expectation value is controlled by a substantially more complicated 
action. Section IV is devoted to these complications and the presentation of the quadratic approximation to the 
action for the current. The case of interacting electrons and photons is considered in Section V. The use of classical 
expectation values, obtained with resolutions belonging to the quantum domain, in particular decoherence, radiational 
time arrow and vacuum polarization are commented on briefly in Section VI. The systematic study of the dependence 
of the expectation values on the space-time resolution around the quantum-classical crossover can be achieved by 
means of the extension of the renormalization group strategy for the density matrix. Some qualitative remarks about 
the renormalized trajectory of this scheme are presented in Section VII. Finally, Section VIII is reserved for the 
summary of the results. The technical parts of the calculation are collected in appendices. The two-point function of 
a hermitian local operator of the CTP formalism is introduced in Appendix A and some details of the calculation of the 
two-point vertex function in the non-interacting Dirac-see are collected in Appendix B. In Appendix C the expressions 
needed for the connected two-point functions for the current and the photon field are presented. The calculation of 
the effective action for the expectation values of these operators is outlined in Appendix D. The evolution equation 
for the functional renormalization group in the CTP formalism is derived in Appendix E. 



II. EXPECTATION VALUES IN QUANTUM FIELD THEORY 



Our goal is to establish relations among the expectation values of a set of local observables like (^'(t)|Oa(a')|^'(i))s 
given in the Schrodinger representation. The problem is interesting when the system is in an excited state. The only 
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technical restriction on the nature of the excitations is that the initial state should be the result of an adiabatic time 
evolution generated from the ground state by adding the source term '^^ja{t,x)Oa{x) to the Hamiltonian density 
where the summation is over the observables considered. Closed systems will be considered below and the issue 
of relaxation and equilibrium in the presence of a reservoir will not be addressed. The state vector of the system 
= U{t,ti)\'i!i)s is given in terms of the initial condition imposed at time tt and the time evolution 

operator U{t,ti) = exp— j(i — ti)H involving the Hamiltonian H. Such expectation values are usually given in the 
Heisenberg representation, (4'j|^(x)|\E',)iif , where x = {t,x) and A{x) = U{ti,t)A{x)U{t,ti) and their perturbation 
series is obtained in terms of Green functions which are different from the usual ones occuring in the expressions of 
scattering amplitudes. Although the trajectory is fixed in the space of states by the initial condition \'^i)s 

the superposition principle allows us to project the system at time on a given final state and to define the 

scalar product {^f\^{t))s = f\U{tf,ti)\^{t))s interpreted as the transition amplitude between initial and final 
states. Its pcrturbative series contain matrix elements like {'i' f\U{tf,t)A{x)U{t,ti)\^i)s which are not intcrpretable 
in terms of measurements according to standard rules of Quantum Mechanics. Such matrix elements are reduced to 
an expectation value for the trivial case = = |0) only (by setting the ground state energy to zero). 

Let us return to our problem, the calculation of the expectation values in excited states and consider the general- 
ization of the forward scattering amplitude for several observables, 

{^i\T[U{tf,U)A„{x„)■■■Al{xl)]\^',)H^{^^\T[Uitf,xl)A„ix„)U{x''„,U)■■■U{U^ (2) 

which differs from the Green function defined as the expectation value of the operator T[An(x„) ■ • ■ Ai{xi)]h, [1], 

{^i\T[Ar,{xn) ■ ■■A,{xi)]\^i)H = {^i\U{ti,tf)T[U{tf,x°JAr,{xn)U{xlti) ■ ■ ■ U{U,x°)A,{x,)U{x'i,ti)]\^i)s. (3) 

The difference occurs in the way the scattering matrix U{tj,ti) appears in these expressions. In Eq. (2) the time 
evolution is constrained in such a manner that the system arrives at the state |\E'j) at time t/, in contrast to Eq. (3) 
where the time evolution is open ended, without any constraint on the evolution after the action of the observables. 
The fimctions given by Eqs. (2) and (3) arc identical for = |0) because the vacuum is stable during the time 
evolution. But the dynamics of expectation values in an excited state, l^'i) |0), requires the use of the expectation 
value (3) obtained in the closed time path (CTP) formalism [2-12], rather than the contribution (2) to the scattering 
amplitude. The careful reader may object that it is sufficient to use the vacuum state because all physical states can 
be obtained from the vacuum by applying local excitations, ie. = Ai{xi)\0)H where Ai{x) is a local operator 

in space-time. Extend the time evolution in this case in such a manner that the system starts at time t = to < ti and 
ends at time t = ti > tf with the vacuum state and write Eq. (3) as 

{0\T[Ai{xi)U{to, h)]T[U{tu x"jAn{xn)U{xZ, ti)--- U{ti,x'i)A,{xi)U{x'i, ti)Ai{xi)U{ti, to)]\0)s, (4) 

where T denotes anti-time ordering. The comparison of this expression with (2) shows that the basic difference 
between the CTP formalism and the usual scattering amplitudes is the presence of the anti-time ordered piece in the 
former case. The proper treatment of such matrix elements requires the introduction of independent time variables 
for the time ordered and the anti-time ordered operators. 

To compare the expectation values obtained by constrained and open ended time evolutions we slightly generalize 
the CTP formalism and allow that the system be described by the density matrices pi and p/ in the initial and the 
final states, respectively. Furthermore, we introduce two sets of external sources, ja{x), coupled linearly to a number 
of local observables, Oa{x), in the hermitian Hamiltonian which is extended to 

H^{t)=HTY. [ d^xjt{t,x)Oa{x). (5) 

Finally, we define the generating functional for the connected Green functions 

^w\j+,j-] ^ *^"W]p/T[e-^4' "''"^^''^pi (6) 

given in the Schrodinger representation. The sources play a double role made explicit by the parameterization 
= j/2 ± j. On the one hand, the variation of j is used to generate the desired expectation values, and on the 
other, j drives the system adiabatically from the vacuum at time t = Iq to the desired initial state at t = ti. 
This allows us to set pi — |0)(0| at to = — oo with the price of having non-vanishing physical external sources, j, for 
some initial times t < ti. The quantity (6) is the transition probability between the states specified by the density 
matrices pi and pf for physical external sources. Expressions like (6) have already been used a number of times with 
Oa chosen to be projection operator, eg. for the time-symmetric formulation of Quantum Mechanics [13, 14] and for 
the study of histories [25-27]. 
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It is natural to introduce the parameterization 



for the field variables. The measured expectation values of the field at a given time t involve the diagonal part of the 
functional {4'~^\pi{t)\4'') with (j)^'^ = 0. The canonical momenta, represented by the operators 11='= = —id/64>^, use a 
slightly extended domain, (f)^'^'^ w 0, of the density matrix. 

There are two time axis in the CTP formalism and each degree of freedom exists in two copies, realized first 
in the case of thermal equilibrium, pi = exp—H/T, [41-45], but can also be clearly seen in the path integral 
representation 

given as a functional integration over the trajectories c/)'^ , cr = ±, with the density matrix elements {(t)'^\p\(t>~) = 
p[(l>^ , 4'~] in the integrand. The duplication is made explicit by defining two kinds of averages, one for each time axis, 

{f[A],T[B]y^ = TVf '^'"~^%fT[Be-'^'^' *''^^^*'^]pi. (9) 



IPS 



Notice that {'1,T[B])^'^ is real according to the relation 



{{T[A],T[B\Y^l^J = {T[B],T[A]);^% (10) 

showing an inherent time inversion invariancc ic. the time inversion leaves Pi, Pf and j invariant and flips the sign of 
the book-keeping variable j only. The time-dependence of the density matrices 

Pf{t) = T[e*-^*''*^W]p/r[e-*-^*''*'^^(*')] (11) 

allows us to write 

eiwU^J-] = Trpf{t)p,it), iU<t< tf). (12) 

The right hand side of this equation is the scalar product of the inital and final hermitean density matrices taken at 

an arbitrary time indicating that the CTP formalism is based on the transition probabilities rather than amplitudes 
and is time reversal invariant. The duplication of the time variables represents the independence of the quantum 
fluctuations in the bras and the kets of the density matrices. We shall argue below that these fluctuations become 
correlated and the usual single time axis formalism is recovered in the macroscopic limit when decoherence [21, 24] 
suppresses the off-diagonal elements of the density matrix. 

It will be instructive to follow the dynamics in the presence of two different boundary conditions in time. The fixed 
boundary condition (FBC) for pure initial and final states, defined by pi = pf = decouples the 

dynamics of the two time axis. The additivity of transition probabilities yields 

=(i,a)°(i) (13) 
p}'+p}',Pi Pf ,Pi Pf ,Pi 

for Tvpf^pf^ = 0. The open boundary condition (OBC) defined by p, = and pf = J2n 1*^"^) (^'("^ | = 1, 

||^("-)^} being a basis for the states corresponds to unconstrained time evolution and couples the dynamics of the 
two time axis at t = tf. Note that the time tf at which the final condition pf = 1 is imposed can be chosen 
arbitrarily for unitary time evolution as long as it is later than the time for which the latest observable inserted. We 
call a boundary condition reflecting if the time at which it is imposed as an initial or final condition influences the 
measurable expectation values. The open boundary condition is non-reflecting because pf = 1 commutes with all 
observables. 

The main virtue of the Heiscnberg representation and the CTP formalism in particular, is to render the initial 
condition problems of Quantum Mechanics simple. The hyperbolic Schrodinger equation allows us to solve the initial 
condition problem in terms of state vectors but these are in general complicated objects. We can deal more efficiently 
with numbers, for instance matrix elements and therefore it is advantageous to convert operator equations into c- 
number equations. But a matrix element contains a bra and a ket, corresponding to the final and initial states. 



5 



respectively as in Eq. (2) and it is not clear how to express the solution of the initial condition problem of the 
Schrodinger equation with open ended time evolution in terms of such matrix elements. This is the problem which is 
solved in the framework of the CTP formalism by extending p, and p/ to positive semidefinite hermitean operators 
beyond the domain of density matrices. 

The direction of the time is determined by the phase of an eigenstate of the Hamiltonian as a function of the time, 
and the time runs in opposite directions along the two time axis. Due to the presence of both time directions there 
is no explicit time arrow in the expectation values which contain both retarded and advanced effects. But one can 
separate these effects, as far as the external sources are concerned, in the case of open boundary condition and we 
shall see below that cf) and (j)'^'^'' of Eq. (7) are the dynamical variables which are build up by retarded and advanced 
effects of the external sources, respectively. 

An important property of the open boundary condition with unitary time evolution is that the generating functional 
is vanishing for physical external sources, W[j^ —j] =0. As a result we have to vary non-physical sources, j, to generate 
the measured expectation values {Oa{x)) = {l,Oa{x)) = {Oa{x),l). This suggests the slight generalization of the 
parameterization of the sources, j ^ j + kj, ie. 

i^ = ^(i±«)±J, (14) 

giving 

^^-^^^^^ = = ^''(^^^ + 1) (15) 

OJa\X) Ij^o ^ ^ 

for arbitrary choice of n. The K-independence of the expectation values reflects a one dimensional symmetry of the 
physical sector of the CTP formalism. 

The parameter k appearing here mixes the non-diagonal quantum fluctuations into the expectation value of the 
observables. We need k 7^ ie. both kind of fluctuations, to derive variational equations of motion for the measured 
expectation values. This condition seems natural because both canonical variables are needed in the dynamics but 
the source j is coupled only to the diagonal quantum fluctuations for k = which are not sufficient to reconstruct 
the expectation values of the canonical momenta 11. The detailed argument goes as follows. The action constructed 
for the field variable </> is obtained in terms of connected Green functions. The n-point functions give rise to 0{(f>'^~^) 
terms in the equation of motion. The simplest linearized equation of motion arises from some connected two-point 
functions. These two-point functions must include field variables at different times to generate time evolution for the 
expectation values. Let us consider the combination 

G = ^(*, I x)[<^(t,y) - <j>{t - At,y)]|*i)H (16) 

of two-point functions with small At as a typical term which can also be written as 

G=-'-{^im,x)m,y),H\\^i)H + 0{At) (17) 

according to Eq. (3). The kinetic energy n^(t, a;)/2 in the Hamiltonian density yields 

G = {^Mt, x)n{t, y)\^,)H + 0{At), (18) 

and shows that the appearance of the canonical momentum operator is unavoidable in the Green functions involving 
fields at different times. We have the option of keeping in the formalism the expectation values of both canonical 
variables, (f> and H, the resulting variational equations being the quantiim analogies of the Hamilton equations. But 
once we have decided to retain the coordinate averages only we need k 7^ to couple both diagonal and off-diagonal 
fiuctuations to the source j. The off-diagonal fluctuations drop out from the field average (15) but the O(j^) terms 
in the generating functional W[j', j] and the 0{(jP) term of the action will retain them. 

A simpler but more formal reasoning starts with the identity j~] = —W[—j~, obtained by comparing 

with its complex conjugate in Eq. (6). This identity reads as j] = — for k = 0, indicating 

that the real part of j] which will be important for the equation of motion is an odd functional of j and has, 
in particular an O(j^) part. As a result the equation of motion for (p contains even powers of the field variables. In 
order to have linear part in the equation of motion we have to allow n ^ 0. 

As a simple demonstration that the measured expectation value results from the retarded effects of the classical 
external sources we consider the linear response formulae. A physical external source, (a;) = —j^{x) = j°'{x) drives 
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the time evolution of the expectation value of the operator Oe{x) in the perturbation expansion of the external source 
according to 

x(I:/^w')^)'^'"'«^->,.=_,-=, (19, 

The linear response formula for Oi{t) is the 0{j) part of the right hand side, 

i f dx'Y^Ux') ({Oa{x'),0^{x))\l^ - {t,T[0a{x')0e{x)])\°J , (20) 

•'a 

which can be written as 

i f f'dt'Y,Ut',x'mt'-t){{Oa{x'),0,{x))\°^^-{^^^^ 
J Jti a 

+i f d^x' f' dt'Y,Ut',x'W-t') ({Oa{x')Mx)fl^^ - {l,T[0^{x')0,{x)])Y^^) . (21) 
J Jti ^ 

If the final time boundary condition is non-reflecting and can be imposed at any time superior to t without modifying 
the expectation values then the relation 

e(i' - t){Oa{x'),Oe{x))l°^^ = e{t' - t){l,T[Oa{x')Oe{x)])ll^ (22) 

holds. In fact, one can reduce tf down to t' followed by the commutation of the perturbation Oa{x') with pf and the 
placing of it at the end of the anti-time ordered part of the expression. This identity shows the cancellation of the 
advanced part of the causal propagator and explains that the usual retarded response formulae are the result of the 
interference between the dynamics evolving along the two different time paths. Whenever a final boundary condition 
is used which can not be displaced in time because it does not commute with all operators in question then advanced 
effects arc left behind. Notice that this remains valid in any order. In fact, the possibility of reducing tf until it 
reaches the time of observation, t, removes any influence of the sources on the expectation value after the time of the 
measurement. The condition of the cancellation was that the non- vanishing sources are classical, j+ = The 
non-classical part of the source, j, representing a coupling among the degrees of freedom of a closed quantum system 
may induce advanced effects. But this component is suppressed in the macroscopic limit by decoherence. 



III. PHOTONS 



The simplest context in which tlw^ dynamics of the field expectation values can be studied is the case of free photons. 
The expectation values are obtainc;d in two steps. First the generating functional for the connected Green functions 
of the conserved current is calculated. This produces the expectation values in terms of the external sources which 
drive the system to the desired initial condition. In the second step a number of different effective actions and their 
variational equations of motion are constructed by performing a Legendre transformation on the external sources. 



A. Photon propagator 



The generating functional for the connected Green functions is constructed by coupling an external current to the 
photon field and using the pair of extended Hamiltonians H — > i?^(t) = J d^xj^'^{t,x)Ai^{x) in Eq. (6). It 
is advantageous to carry out the calculations in the path integral formalism where the action for the photon field is 
written as 



Sm[A] 



L 



(23) 



7 



Here Do is the free photon propagator Dq ^ = Dj.^ + D^^ and the projection operators onto the transverse and 
longitudinal components of the photon field T^" = g'^'^ — L'^" and L'^'^ = d^^d" /H are used to construct D:^^ = (□— ie)T 
and Dl^ =^{n-ie)L. 

A few words about the notations: The field configurations (p{x) are usually considered as vectors, (jjx, and space-time 
integrals with occasional summation over repeated indices as scalar products, eg. / dx(j){x)x{x) = J 4>xXx = </> • X 
or J dxA^{x)j^{x) = A - j. The space-time coordinates are written as a; = {t,x), the space-time indices are given 
by Greek letters, the Latin letters denote combined indices like a = {n,x). Repeated indices are summed/integrated 
over. 

The generating functional is given in the framework of the path integral representation as 



j D[A]e^'^-^o'-A+^3-A^ (24) 



where the two component integral variables A = A~) and external sources j = [j^ are introduced to simplify 

the expressions. The current will be parameterized as in Eq. (14). Some complications arise due to the presence 
of the density matrices pi and p/ in the functional integral. Anticipating the perturbation expansion we use Gaussian 
density matrices which generate quadratic expressions in the action of the path integral, ie. In/>j[j4_|_^t^; and 
\n.pf[A-^tf\Aj^^tj\ should be at most quadratic in the fields. The inverse block propagator is 

v=(Y ^^^^ 

where the second term on the right hand sides stands for the contributions of the density matrices and is non- 
vanishing for t = U OT tf only. Notice the presence of an infinitesimal imaginary part in the inverse propagator 
(25). The operator expression (6) contains a hermitian Hamiltonian because the (anti)time ordering implies the time 
boundary conditions instruction. In contrast, the time boundary conditions involve the infinitesimal imaginary parts 
in the inverse propagators in the path integral representation. 

Most of the work reported below is based on OBC with the vacuum as initial state. Although one can determine Dq^ 
directly in the path integral formalism for this boundary condition it is simpler to construct the photon propagator, 
Do, within the operator formalism. The result is 

° {^Dl^' + i^Do -D'^'^'^' + i^DoJ ^ ' 

in the base (j^, j~) where and D*^" denote the near and far field Green fimctions [46], respectively and i^Do 

stands for the imaginary part of the causal (Feynman) Green function. The retarded and advanced Green functions, 
£)ret ^ ]jnei.v ^ i£,far ^^^^ £,adv ^ £)near _ i£,far frequently be used, as well, cf. Appendix A. 

The Gaussian integral can easily be carried out in Eq. (24) yielding 

W^''"'\j] = ~j-Do-j (27) 

and a straightforward calculation gives 

W^^'^'lil]--^-''^ .( ^ ^o'^ ].(A (28) 

i^OBcLJ'JJ- 2 V£>5<'t kD'^''''' + i^DoJ \jj ^ ' 

for open boundary condition. In the case of fixed boundary condition the two time axis decouple yielding 
W^^Bcbl = -lj+-{Dr' + i'^Do)-j+ + lj--{Dr'-i'^Do)-j- 



±1 " 

2 



2 \D'^<"'' + Ki'^Do kD^""^' + ^i'^DoJ \j 



(29) 



B. Effective action 



The role of the effective action is to provide the functional for the variational equations satisfied by the space- 
time dependent expectation values and it is introduced by performing a functional Legendre transformation on the 
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generating functional of connected Green functions. Special attention must be payed to the fact that the Legendre pair 
of the sources, the expectation values, are complex in general. Accordingly we perform the Legendre transformation 
separately for the real and imaginary part of W[j] = ^W[j]+i?sW[j]. The measured expectation value of the photon 
field is always real and it will be determined by Therefore, we start with the eflfective action 

pphot ^adv] ^ SRI^phot J] _ J . ^adv _ j . A, (30) 

with independent variables consisting of the measured expectation value of the photon field 



and an auxiliary field variable 



(31) 



^adv ^ (32) 



The inverse Legendre transform based on the relations (30) 



' = - SA 

jrphot[^^^advJ ^^^^ 



serves as equations of motion for the expectation values. Since the averages A and A'^^ receive contributions from 
the diagonal and non-diagonal fluctuations, respectively, these equations of motion control the time dependence for 

both the coordinates and the momenta. 
The generating functional (28) yields 

T^^cl^, A"'^^] = -A'"^^ ■ • A + ^A^^^ • • D'^""' ■ D^^''^ ■ A'"^^, (34) 

for open boundary condition. The corresponding equations of motion for A and A^^^ are 

^adv ^ £,adv^.^ (35^ 

and 



A = Dl^'j + KDl'^^'^j, (36) 

respectively. They show that non-diagonal fluctuations contribute to the advanced fleld A^'^ generated by j and the 
physical expectation value. A, is indeed the retarded field created by the physical external current j. The two fields, 
A and A^'' remain decoupled in the absence of interactions. For fixed boundary conditions, Eq. (29), we have 

pphot yl^"^^] = ^adv _ jyaea.r—1 _ ^ _|_ ^^adv _ ^near— 1 _ ^adv (37) 

This expression can be obtained from Eq. (34) by the replacements D}y^ — > Dq^^^ and Dq^'^ — > Dq®'^'' which represent 
the loss of the interference between the two time axis. 

The eff'ective actions introduced so far provide the equations of motion for the expectation values for both field 
variables of the CTP formalism. But the field A^"^^ is not physical and its presence is not necessary to extract the 
time dependence of the photon field A. We simplify the Legendre transformation by keeping j as a fixed parameter. 
The resulting effective action for the photon field alone is 

rP^^°^ [A] = 3?W^P''°* [j, ~j] -j-A, (38) 

where the dependence of rP'^°*[A] on j is not shown explicitely. The effective action for OBC turns out to be 

= -^(^ - i^o'^) • ^0^'"""' • - ^ri), (39) 
generating an equation of motion (5rP'^°*[yl]/5A = — j identical to Eq. (36). 
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So far the value of the parameter k entering by expression (14) is arbitrary. The choice k = seems to be natural 
and simple and the corresponding currents j and j generate retarded or advanced Lienard-Wiechert potentials, 
respectively, according to the equations of motion. But the effective action (34) contains no 0(j4"^) piece when k = 
in agreement with the general remark made in Section II stating that the non-physical field can not be eliminated 
and no classical action and variation principle can be found for the measured expectation value of the photon field. 
This is reflected in the appearance of the coefficient 1/k in the effective action (38) because for k = the generating 
functional (28) has no 0{j^) term in the real part and the Legendre transformation is not defined for linear functions. 

For OBC the Legendre transform of the imaginary part is defined as 

r&BC '"[^™] = S^OBC \J] - J ■ (40) 

where 



The corresponding equation of motion is 



3 = "Xl^J • (42) 

The generating functional (28) yields the explicit form 

roBC ™ [^™] = - • ^^0 ^ • (43) 
For FBC the functional depends on both currents and the Legendre transformation results in 

rP^°^ ™ [A™ , A""^" ™] = 9 W^Ibc j] - J • ^™ - J • ^""^^ ™ (44) 

where 



^adv im _ <^^^FBcb'i] 



^adv im _ tuij J i (45) 

^3 



The actual form, 



, . 1 I ^2 

rP^°^ ™[A™, yl^dv im] ^ I^^adv imc-^J-l^adv im _ ^!mcj£,-l^im _^ ^^imcj^,^ l^adv im^ ^^g^ 

follows from the functional (29). 

These effective actions will be used in Section VIA. 



IV. CURRENT DYNAMICS IN THE NON-INTERACTING DIRAC-SEE 



For a non-interacting system the determination of the expectation value of a local observable is trivial as long 
as the observable is a one-body operator, as in the case of the photon field discussed in the previous Section. The 
triviality comes from the fact that the generating functional for the connected Green function of the elementary fields 
is quadratic and can exactly be calculated. Higher order Green functions, {0\T[4ixi ■ ■ ■ 0x2,i]|O) with n > 2, factorize 
to a sum of n disconnected products of two-point functions, according to Wick-theorem. The physical origin of this 
factorization is obvious, it comes from the absence of interactions between the particles. Each particle created by 
one of the operators must be destroyed by another operator in order to find a non-vanishing contribution to the 
vacuum expectation value. In other words, the largest cluster with non-factorizable structure contains two operators. 
But we meet serious difficulties as soon as the expectation value of such an operator is sought which controls more 
than one particle. The non-factorization arises from the possibility of creating particles by one operator which are 
annihilated by different other operators. One has to perform a "bosonisation" in a fermionic system because only 
bosonic operators can have non-vanishing expectation values. The simplest and most important bosonic operator is 
the electric current j considered in this Section. 
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Let us first introduce the generating functional, T^^'^'fa], for the connected Green functions of the electric current 

by means of the functional 



where the two-component fields -ip = {tp+jtp-) a = (a+,a ) were introduced together with the inverse electron 
propagator 

^o'=(^S _^OG-it^o)+GBC> (48) 

where Gq^ =ip—m. We shall need later the composite operator counterterm 

AZ3 + /3^ [D-^ 



ocrm = « • ) ■ « (49) 

with AZ3 being UV divergent and the finite part /3 being fixed by a renormalization condition. 

In the absence of charges in the initial and final states we identify the generating functional as W^o'[a] = M^'^'[a, 0, 0]. 
When the system contains n~ electrons and n+ positrons we take pi = or pf = 

I*) = fl ( / 11 ( / xL^*..) 10), (50) 

with t = ti or tf, respectively. The wave functions x~ and x+ describe the one-particle states. The generating 
functional for the current is then written as 

We want to generate the expectation value of the current 

i^? = ^fe7''^.-(7'*^.)*'"^r] (52) 
which changes sign under charge conjugation. To this end we use the 7-matrices 

where = 0+ in the Dirac Lagrangian and the minimal coupling. In this notation the current is ja = '4> ' la ' V'- 
The mathematical source of the complications is that 

W^a] = -iTr In G-i[o] + SctIo] (54) 

where (^^^[a] = G^^ + ^ is an involved functional even for the vacuum, ie. in the absence of initial or final charges. 
It is easy to understand the origin of this complication in the special case of pi = p/ = |0)(0| and time- independent 
source a. The generating functional T^'^'fa] is now the difference of the energy of the Dirac-see in the presence of 
the vector potential a"*" and a~ . A local modification of the vector potentials creates a non-local polarization of 
the Dirac-see with all filled one-electron states contributing. Therefore, the infinitely many negative energy states 
filled up in the Dirac-see renders the functional highly non-trivial. It is not necessary that the states be filled. For 
scalar particles the structure of the generating functional remains the same except of the change of the overall sign. 
The external potential coupled to the particles creates a polarization of the ground state which involves an arbitrary 
number of (non-interacting) particles. The expectation value of bilinear operators measures the polarization created 
by the given boundary condition in time and its space-time dependence satisfies highly non-trivial equations reflecting 
the multi-particle dynamics of polarization in the ground state. 

Our strategy followed in the case of the non-interacting Dirac-see will be similar to that of Section III, ie. we 
calculate first the generating functional for the connected Green functions of the current and then proceed with the 
construction of the effective actions. For weak fields the functional can easily be calculated in the framework of the 
perturbation expansion, which is the main task of this Section. The issue of strong fields is far more difficult and will 
only be commented on briefly. 
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A. Connected Green functions 



We start with the simpler case when there are no additional charges immersed in the Dirac-see. For weak external 
sources the functional (54) can be written as a functional Taylor series, 

OO -J 

^*^'[«] = E;i! (55) 

Tl=l 

where the super-index a = (±, x, /x) identifies a time axis, a space-time location and a vector index, and the coefficients 
give the connected Green functions. The expansion of the logarithmic function in Eq. (54) results in 



W''[a] = -iTi In Go'+iJ2 H7^Tr(Go • ^)" + Scria] (56) 

n=l " 



and 



where 



= -^zr- E Tr[Go • 7l,(,, • Go • Tip,,, • • • Go • 7t(„) • Go • 7|] - Sn^^Z, + P)Dt (57) 
Pes„ 



r. = [l- . (58) 

The symmetrization with respect to the exchange of the external legs is achieved by the summation over the permu- 
tations of the vertices in (57). The odd orders are vanishing according to Furry's theorem. For the sake of simplicity 
we truncate the generating functional at the quadratic order and write 

W^\a] = ~a-dR-a (59) 

where the renormalized current two-point function Gb = Go — PDt is given by 

Gj+6 = -^Tr[G++ . 7„ . G++ • 7^,] + AZsT{a - ie) 
Goffe = -^Tr[G+-.7a-Go-+-76] 

Goa'b = -*Tr[Go"" • 7a • Go~ • 76] - AZsT{a + ze). (60) 

The notation G = G"*""*" will be used in the rest of the paper. Well known results (eg. Ref. [47]) include the 
renormalized two point function 



= / ' 



G^; = / e-''-G^:'o 



157r 



(61) 



obtained in the framework of the gradient expansion where SJG = since the creation of a mass-shell e e+ pair is 
forbidden. We shall use the same parameterization of the propagator 

~ ^ far'' + i^Go -iGf- + i^Go] , . 

° y^Gi^'^ + i^Go -G^'^'^' + i^GoJ' ^ ' 

as for photons and the retarded and advanced current Green functions will be defined by Gq''* = Gq^^^ + \G^o^ and 

Gf^ = G;J'=^'' - \G\f., respectively. 

An interesting feature of Eq. (57) is the need of renormalization. The Green functions VFjf'aj are finite for 
n > 3 but the two-point function diverges. In fact, this two-point function is identical to the one-loop photon self 

energy except of the missing factor e^. The electrons of the Dirac-see tend to approach each others too frequently 
and make the two point function divergent when the two legs approach each other in space-time. This is a well 
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known problem in QED and the cure is the introduction of the counterterm SctIo]- The UV finite part. /?, of the 
counterterm is fixed by a renormahzation condition to be imposed. The lesson of this divergence, a relativistic effect, 
is that the dynamics of the current j can not be defined without an additional scale, the cutoff, even in the absence 
of interactions. 

The need of a renormahzation condition for non-interacting particles demonstrates a characteristic difference be- 
tween first and second quantized systems. In quantum mechanics observables are defined by the operators. In quantum 
field theory however the observables represented by composite operators may need counterterms and their proper def- 
inition must include the corresponding renormahzation condition. Symmetry principles can not fix the counterterm 
to the current as long as it is transverse and we find a one-parameter family of current with the two-point function 

^(0T[J•^J-,]0|)+/3□^,,,,T''^ (63) 

as far as the many-body aspects are concerned. The finite part of the counterterm influences the product of two 

current operators at zero separation and this contact term plays an important role even at finite energies. We shall 
see later that the current coupled to the photon field is defined by /? = 0. 

The generating functional is more complicated in the presence of additional charges in the Dirac-see. We then have 
the additional term ■ G[a]o ■ fj in W[a, fj, fj] which gives for instance 

W^'la] = W^|Bc[a] - i E HGl^]uf,.,),Ui^..)) (64) 

<T = ± 

for charged fixed boundary condition with a single electron. The initial (final) state is characterized by the space-time 
point Xi {xf) with bispinor index ji {jf) and a = {x,j) denotes the combined index. As an other example, the open 
boundary condition for two electrons with a„, n = 1, 2 in the initial state yields 

Now the generating functional contains odd orders in a. The different structure of the see and valence contributions, 
the first and the second terms on the right hand side of Eqs. (64)- (65), reflects the fact that the Dirac-see is made up 
of negative energy one-particle states while the additional charges, introduced by the creation operators, correspond 
to positive energy. As a result, the valence charges move freely while the motion of the particles making up the 
Dirac-see is restricted by the Pauli-blo eking. 



B. Effective actions 



The equation of motion for the current expectation values is derived from the effective action 

pel J jadv] ^ ^el [a,d]-d- J'"^^ -a- J (66) 

involving the sources 

= ± a (67) 

and independent variables 



J 



5a 



J-v = ^l^phA (68) 
oa 

which are conserved currents. 

The effective action for the physical fleld A only is defined by 

T''\j]=mV''\a,d]-a-J. (69) 

The actual expressions for the Dirac-see can be obtained by the replacements Do ^ G, A — !■ J, and A^'' J*^^ in 
Eqs. (34), (37) and (39). 

When charges are added to the Dirac-see the valence propagator contributes to the a-dependence in VF®^[a]. This 
requires j ^ and a redefinition of G in Eq. (59) to be used in the Legendre transformation. 
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C. Polarized charges 

Let us consider now the dynamics of localized charges polarized from the Dirac-see and described by the effective 
action truncated at the quadratic order, 

= -lj.r(2>'- J. (70) 

The current is supposed to be slowly varying over distances £ < 1/m rendering the leading order gradient expansion 
(61) 

r(2)el = T. (-71) 

applicable. Feynman's e-prescription is displayed in the expression explicitely for the calculation of the inverse of the 

kernel. This expression obtained for OBC holds for FBC, as well, as long as the total charge is vanishing and 
consists of closed flux-tubes. The calulation indicated briefly in Appendix B for the choice /3 = 0, justified for the 
electric current in Section V A, gives 

j.(2)el..^_WQ^^^_^,^2^y,.^ (72) 

ie. the vertex function is step function, assuming the value or — 15m^/8 for spacelike or timelike separations, 
respectively. 

Let us now assume the form 

J^f =5''°^7V,p(|x-x,|) (73) 

3 

for the current where Nj denotes the number of electrons of a localized state and p{r) = 1. Notice that the charge 
is polarized out from the negative energy one-particle states for weak external source and there is no reason for Nj 
to be integer. The average distance between two world-tubes is given by 

rjk= p{\y-Xj\)\y- z\p{\z-Xk\). (74) 
To identify the total mechanical energy we need E'tot = obc['^]/(*/ ~ *») ™ limit tf — U ^ oo, 



■E'tot = 



15m 



-{tf-ti)Nl^ + Y,NjNkrjk 



(75) 



A'tot = ^3 being the total charge. The effective action for /3 = is IR finite for neutral systems only. This is a 
trivial result, reflecting the impossibility of polarizing out a net charge by a neutral external source. 

It is interesting to note that the action of the current defined by fi < is that of a classical electrodynamical action 
with charge g = 1/V— /3 and electric susceptibility 0{g^k^ /m^), cf. Section VIC yielding 

^t°t = -|-E^i^M Pi\^-^3\) r— n P{\y-Xk\). (76) 

°^ Jx,y y\ 

The screening efii'ect of vacuum polarization detected by the current defined with /? < removes the IR divergence 
and allows a net charge. 

Though it is perplexing to find interactions between localized charges in the non-interacting Dirac-see there is a 
simple explanation for it. In fact, let us consider the normalized two-particle state 

ll){xi,X2) = —^[lpi{xi)ll)2{X2) - '4'l{X2)lp2{xi)] (77) 

in non-relativistic Quantum Mechanics where tpi{x) and ipi{x) are two not necessarily orthogonal wave functions. 
The matrix element 

{^\A\^) = i^[(V'i|A|Vi)(V2|V'2) + (V2|A|V;2)(V'i|V'i) - 2i?e(V'i|A|V'2)(V'2|V'i)] (78) 
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shows that the kinetic energy receives a contribution from the entanglement of the two particle states in the region 
where the one-particle wave functions overlap and this additional piece in the energy might be interpreted as an 
interaction potential 

(Vi|^i)(^2|V2)-KVi|^2)p- ^''^ 

All this is the well known Pauli-blocking: the anti-symmetrization excludes certain kinematical regions leading thereby 
to an increase of the kinetic energy. This blocking is restricted to small separations since the one-particle states must 
overlap. 

But this picture is static, ignoring the dynamics. The modification of the density distribution can be viewed as a 
result of polarizations due to particle-anti particle pairs whose propagation establishes non-local exchange correlation 
effects during some time. The near-field vertex function obtained from the 0(n^^) kernel of the non-physical current 
with /3 < restricts correlations onto the light-cone. This is what happens in QED, namely the mixing of the states 
of the photon and the electron-positron pair in a gauge and Lorentz invariant dynamics keeps the polarization of 
the Dirac-sec on the mass shell m? = 0. But the the UV regime of the 0(0^^) kernel of the true electric current is 
suppressed and the singularity at the light-cone, ie. at vanishing invariant length is reduced to a finite discontinuity 
leaving no asymptotic particle-like states contributing to the current two-point function. It is remarkable that the 
vertex function is actually a step function and it establishes a distance independent force between static world-tubes 
in sufficiently long time. 

Notice furthermore that it is not necessary to have relativistic fermions to find such a long-range exchange correla- 
tion. Charged relativistic bosons display a similar dynamics, the only difference compared to the fermion case being 
an overall sign in the generating functional W[a] apart of the modification of the current operator the external field 
6 is coupled to. But both the bosonic and the fcrmionic currents arc bilinears of the elementary fields and, therefore, 
the vacuum polarization should be qualitatively similar. It is not the filled states of the Dirac-see which are important 
but rather the propagation of states created by the current operator from the vacuum. These excitations are always 
bosonic and thus subject to an effective bosonic description. Non-relativistic systems, such as non-interacting fermions 
at finite density display similar effect, too. 

This distance-independent force would not be observable even if one had sufficient space resolution in the experi- 
mental device. The linear potential is a special feature of the non-interacting Dirac-see only and becomes screened 
by the electromagnetic interactions. In fact, suppose we had a matter and an anti-matter localized states and we 
started to separate the two components. Increasing the separation beyond the Compton wavelength of the electron 
costs energy comparable with those needed to put a virtual e~e+ pair on the mass-shell. Furthermore, the gradient 
expansion is valid for separations larger than the Compton wavelength. Therefore the charge of the real particles 
which are created by the minimal coupling vertex tends to screen the original charges and only localized states with 
non-integer charge remain strongly correlated. 



D. Localized states 



The preceding discussion of the dynamics of the polarized current is restricted to weak external sources. The 
involved, non-polynomial character of the generating functional VFo'[a] makes its appearance for strong external 
sources. Though the correspondence o ^ J is unique but its inverse is not and the effective actions are not uniquely 
defined for strong external sources. By retaining the profile of the current only one looses important information. For 
instance, a given charge density has different dynamics depending on the nature of the states which contribute to the 
charge. This is reminiscent of the small and large polaron problem in solid state physics where the same polarization 
may induce different dynamical responses depending on the way the polarization builds up. 

For weak external sources the polarized charge is made up by small contributions received form a number of negative 
energy one-particle states. The resulting charge density can be localized but extended states contribute only. As a 
result, Pauli-blocking renders the dynamics of such localized states highly correlated. For slightly stronger localized 
external field a bound states are formed in the mass gap. This case will be considered in a qualitative manner below. 
For stronger localized external field the localized states may appear in the negative energy continuum. Such bound 
states appear as a violation of the convexity of the effective action; a rearrangement of the vacuum takes place. One 
member of a virtual e+e~ pair created by the external field jumps into the potential well and the resulting energy 
is sufficient to put the other charge on mass-shell. As a result, a new, well defined Legendrc transform of VF(f [a] is 
regained. In other words, the effective action, r'''[J], is multi- valued and may have several consistent "sheets" for a 
given current. 

Let us return to the case of medium strong, localized external field which create localized states in the mass gap. 
The external fields fi and fj can be used to place charges into the initial or final state which can be captured by these 
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potential centers. These charges have weak overlap with the filled extended states and are supposed to obey the 
"free dynamics" anticipated from classical mechanics. We start with charges of the same sign with fixed boundary 
condition where W'"'''[fl] = M^[a]3 — ilnG'[a], G[a] being the valence propagator with given initial and final points. 
We are interested in the possible decoupling of the trajectories. Therefore, the initial and final location and spin 
of the charges are chosen to be the same. The external potential is chosen to be static and of the form = g^ow 
where the temporal component, u, is the sum of potential wells which are strong enough to create localized state in 
the mass gap and spread enough such that the characteristic size of the bound state, £, be large compared to the 
Compton wavelength of electrons, mi > 1, in order to avoid pair creation. We choose tf —ti^ 1/m and shall apply 
the non-relativistic approximation for the valence propagator. What is crucial is that T4^[a]o can then be ignored in 
W^^ [d] beside of the valence propagator. Once the vacuum polarization is suppressed, the rest is a simple problem in 
non-relativistic quantum mechanics. In fact, we have at this point 



for t/ — — > c» where Eq denotes the ground state energy of electrons bound to the potential wells. When separating 
the wells the charge density remains localized at the wells and the contribution of the last term on the right hand 
side breaks off into a sum of the contributions of separate wells. The correlation energy of Eq. (79) becomes small 
for well separated states and therefore the contributions of the potential wells to Eq decouple. 

Once the separation of localized charges is established wc can consider the issue of the dynamics of a single localized 
state. Let us suppose that the current J is a narrow flux tube with flux N and its center follows the world- line x^{s) 
in space-time. One interprets this as the motion of a charge Ne along the world line x'^{s). If the value of the effective 
action evaluated for this current, divided by the invariant length of the world line x'^(s) is independent of the time 
of evolution then this ratio is the mechanical mass. Let us consider a world line without acceleration. Then Lorentz 
invariance allows us to set Xi = xj. We use the same external potential as in the previous case except that it has a 
single potential well only and we place a single electron in the initial state. The generating functional for this problem 
is (64) whose value can be written as in Eq. (80). The ground state energy, Eq, is the sum of the expectation values 
of the kinetic and the potential energies. The latter cancels against the last term on the right hand side of Eq. (80) 
because J^{t,x) is the particle density, ie. the magnitude square of the ground state wave function apart of some 
finite time at the beginning and at the end of the motion which give negligible contribution in the limit U = tf ^ oo. 
We are finally left with r'''[J] = — {tf — ti){m + Ek) and the mechanical mass, the coefficient in front of the first term 
on the right hand side of Eq. (1), can be written as 



where Aruk is the expectation value of the kinetic energy in the bound state which is approximately l/(£^m). The 
second term represents the effects of the polarizations of the Dirac-see by the external potential coming from M^[a]g. 
Its order of magnitude, m£u, is estimated as the product of the polarized charge square, assumed to be proportional 
to u/m, the ratio of the energy scale of the potential and the mass, the length scale of the potential well £ and the 
string tension a w introduced in Section IV C. Finally, Arrij stands for the radiative corrections arising from the 
electromagnetic interaction and is e^/£. 

The mechanical mass depends on the details the elementary particles localized in their world-tubes. Once the 
coupling to the photon field is introduced the point-like charges develop their vacuum polarization cloud, i « l/m, 
and the correction to the mechanical mass becomes non-perturbative, 0{m), but uniquely defined and can be removed 
by renormalization of the mass. 



After the separate discussion of the dynamics of the expectation value of the photon field and the electron current 
we embark the case of the interacting system. 




(80) 



M = m + Am/s -|- Ampoi + Am- 



'7) 



(81) 



V. INTERACTING ELECTRONS AND PHOTONS 



A. Generating functional for connected Green functions 



The generating functional for the connected Green functions of the current and the photon field is 




(82) 
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where the minimal coupling contains the matrix 

The pcrturbativc rcnormalization is carried out by using the counterterms Sct = S'^t^ + S'^'rp where S^^^ contains 
the usual counterterm needed for the rcnormalization of the electron and photon Green functions in a given order of 
the loop-expansion. The additional term is chosen to be 

= -{AZ3 - a)ea ■ Dq^ ■ A+ ^(AZa + 0)a ■ Dq^ ■ a (84) 

The origin of this expression becomes clear by considering the photon part of the action, the sum of S^j,, the bare 
action and the photon self energy counterterm, 



Sct + St^ ■Do^-^=\a- D^^ ■ A + ^(ei - a) ■ D^^ ■ (ei - a) + aa ■ D^^ ■ A+ ■ D^^ ■ a, (85) 
by using 



i- = l+AZ3. (86) 

The second term on the right hand side with a = is the counterterm for the photon self energy but its dynamical 

origin is in the fermionic sector of the theory. Any vector potential which couples to the electric current as the photon 
field must appear in this counterterm on equal footing with A. This shows that the choice a = /? = of the finite 
part of the counterterms is needed to identify the current coupled to the source a with the electromagnetic current. 
The electron field can be integrated out easily in Eq. (82), leaving behind 



where [a] stands for the generator functional of the connected Green functions of the current in the non- interacting 
Dirac-see as given in eqs. (54), (64) or (65) for different boundary conditions. The integration over the photon field 

can be carried out in the loop-expansion, see Appendix C for the details. We record here the results for the quadratic 
part of the generating functional in the absence of charges in the initial state, 

where 



G=— (89) 

Go 

is the full current two-point function containing the current self-energy E**^ and 

D = -— -^-^ (90) 

stands for the photon propagator involving the photon self-energy Y,p^°^. The formal expressions for the self-energies 

up to two-loops are 

= e'^aDoa - e^cC ' "^"d^^ ' (91) 

and 

SPhot = e2a(6o - e^W'§''^)a (92) 
where the first and second terms on the right hand sides represent the one- and two-loop corrections with 

<it^=<l'i('^^^oa)e. (93) 
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being the one-loop self-energy and vertex correction to the current two-point function. 
The more useful parameterizations (14), (67) yields the form 



^W[a,a,j,j] = ^{a,j,a,j)- ( '^''^j^^^'^ ^1 • 



(94) 



with 



and 



-K^OBC = — 



-K'fbc 



^near 
near ^ ^near 



e(G 



near y-jnear 
• ^0 



• 9£»o) 



Do 



(95) 



(96) 



in the absence of charges with positive energy, for open and fixed boundary conditions, respectively. As pointed out 
after Eq. (36) the real part of the non-interacting generating functional for fixed bomidary condition can be obtained 
from that of the open boundary condition by replacing the retarded and advanced Green functions by the near field 
version. Such a change is not sufficient for interacting system where higher loop may bring in the product of two 
imaginary parts into 3?W^. However by restricting our discussion to distance scales beyond the Compton wavelength 
of the electron the creations of mass-shell e^e~ pairs are suppressed and 3G = 0, cf. Eq. (61) obtained in the 
framework of the gradient expansion. It is pointed out in Appendix D 2 that this is the only violation of the simple 
rule mentioned above and once the pair creation is excluded one can again obtain easily the fixed boundary condition 
expressions from their open boundary counterparts. 

The comparison of Eqs. (C20) and (C21) reveals that the imaginary part of the propagator undergoes more 
substantial changes then the real part when the boundary condition is modified. 



B. Legendre transformation 



We introduce a number of effective actions. The simplest effective action arises from the Legendre transformation 
of the variables and j^, 

T[J,A\ = W[a,j]-d-J-j-A (97) 

where 
and 

A^'-Z. (99) 

53 

The inverse of the block-matrix needed for the change of variables in the generating functional (88) can be read off 
from the result of Eq. (D3) yielding 

r[J, A] = T"'<'^^[J] + ^A-Do^ -A-eAa-J (100) 
where the first term is responsible for the dynamics of charges moving within the interacting Dirac-see and is given 

by 

on the two-loop level. The equations of motions are 

ffii^ = -a (102) 
SJ 
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and 



6T[J,A] 
SA 



-J- 



(103) 



What is remarkable in the expression (100) is that the electromagnetic interactions are represented without loop- 
corrections, ie. by the inverse of the free photon propagator and the minimal coupling without a form factor. This 
is bccaiisc the introduction of separate variables for the current and the photon field formally removes the one-loop 
contributions from the action. In fact, the elimination of one of the fields by its equation of motion reintroduces the 
one-loop pieces in form of an action at a distance. In fact, eliminating J or A by their equations of motion generates 
the action at a distance effective actions 



r[i] = 



or 



-J-G -J 



(104) 



(105) 



for J = 0, = 0, involving the full propagators. 

Although the effective action (100) shows obvious similarities with the action of classical electrodynamics it is not 
satisfactory. In fact, though either fields of the type + or — can be used to read off measured expectation values but 
we have both of them in the effective action coupled to each other. The physical expectation values are easier to read 
off when the Legendre transformation is performed on the variables introduced by Eqs. (14) and (67). For this end 
we shall use the effective action 



r[J, J'^'', A, A'^^] = ^W[a, a, j, j]-J-a- J^^^ -a-A-j- A'^'' ■ j 



with the new variables 



J 



ad 



V SdiW[a,j] 



6a Sa 
The inverse Legendre transformation is based on the relations 



, _ 5diW[a,j] 



and 



J = -- 



SJ 



5T[J, J'"^'' , A, A""^"] 
SA 



ST[J,J^'',A,A'"^''] 
5T[J, J^'^'' , A, A"^^] 



(106) 



(107) 



(108) 



(109) 



The block-matrix matrix of Eq. (95) yields, after some straightforward calculation outlined in Appendix D 3, the 
effective action 



r[J, J^'^", A, A'"^''] = - 2 (•^' ^' J^'^^^ ^'"^^) 



tr-l 



K 



/ J \ 

A 

jadv 
y^advy 



(110) 



where 



(111) 



for open boundary condition. The retarded part W'^'"'^'^^'^* = ^W^^'^^^'^ + ^W'^'"^'' is defined in the same manner 
as for the free block- propagator in Eq. (26), cf. Appendix D2. The equations of motion (D23) give J*^"^ = ^4*"^"^ = 
for vanishing off-shell sources, a = j = 0, as expected, and we find the equations of motion 



J 



1 + e^W' 



1 

!l(4)r ^ret 



D 



-.Gl^^.{eA-d), 







A = .(eJ-J) 



(112) 
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for the physical fields. The one-loop contributions are again absent and the two-loop term represents the electro- 
magnetic interaction in the Dirac-see in the first line of Eqs. (112) as in the effective action (100). According to the 
second line of Eqs. (112) both the external and the dynamical currents, j and J induce retarded potentials. 

In the case of fixed boundary condition we replace the retarded and advanced propagators by the near field version 
in the absence of creation of mass shell e+e~ pairs and use 

/rYnear-l , „2/^near-l TTrel(4)near /=,near-l „ \ 

K-' = -n '.^^ ~i (113) 



in Eq. (110). The equations of motion become 



G^^'^'a = -(l + e^V^;,'^^)""*''.^^;^^'-!). J + eGS^^"-- A 

D^''^'] = eDS"'"'' J - ^. (114) 
The simplest effective action includes physical fields only and it is defined as 

r[J, A] = ^W[a, a,j,j]-J-a-A-j (115) 

with the new variables 

Sa Sj 

and the source a and j treated as passive parameters in the Legendre transformation. Notice that the important block 
in the generating functional (94) we have to invert to find the equations of motion for the physical field, the upper 
left one, is the same for both boundary conditions. The dependence on the boundary condition appears through the 
linear pieces of the effective action only. The simple steps shown in Appendix D 4 lead to the effective action 

kTIJ, A] = r"''"'^[J] + ^A ■ D'^''^'-'^ -A-eA-J-A- 

pmechjjj ^ lj.(^near-l_^^near-l.^el(4)near_^near-l). j_ j.^ext (n^) 



where the external sources 



jext ^ D'^'"''-'^ ■ Wj - eWa 

^ext ^ (^near-l^g,near-l.^cl(4)ncar^^„car-l).^^_g^^^ (^^3) 

are given in terms of the boundary condition dependent terms Wa, Wj shown in Eqs. (D25) and (D26). The 
corresponding equations of motion for the field J is 

J = -G"<'°^^(a + eOf^' ■ j) + z z [A + D'^^^' ■ j + eOf^' ■ G'^^^' ■ a)] (119) 

^ " -"^ (^near-l _|_ ^near-l ^^(4)n ^ ^near-l ^ J ' U yj V / 

where the upper (lower) indices correspond to open (fixed) boundary condition. The equation of motion for the 
electromagnetic field is 



A = e£»j;°" • J - Dl^^ ■ j + -D'^' ■ G^°* • (a + eL»J,'=' -3) + {e'^Dl^^ ■ W'^'^"^'"^ • D5<=* - e^DJ,'=' • G[,°' • D'^''' ■ G[,°' • Dl^') ■ j (120) 
for open boundary condition and 

A = e£>5^'*' • J - [D"^^"' - e^Dl""' ■ . Dl"^' + e^Dl""' ■ Gg^"*' • Dl"^' ■ Gl""' ■ Dl""'] ■ j (121) 

for closed boundary condition. 

Notice the difference between Eqs. (112) and (119)-(120), the equations of motion for the physical fields derived 
from the actions T[J, J^^, A, A'^'^"] and T[J, A] with open boundary conditions. We shall come back to this apparent 
paradox of having different equations for the same quantity in Section VII. 
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VI. MICROSCOPICAL CLASSICAL FIELD THEORY 



It is important to distinguish two length scales when constructing classical field theories from quantum theory. 
The shorter one, a = 27rA~^, is the UV cutoff of the underlying Quantum Field Theory. The other is the quantum- 
classical crossover length scale, ^crj which is a rough order of magnitude estimate since the scale of the actual crossover 
depends on the environment. Classical field theories derived above for the expectation values of local operators provide 
information about the dynamics with a space-time resolution limited by a only. Structures seen at scales a < ^ < 
characterize the microscopic quantum dynamics. Our classical field theory can be viewed as a time dependent 
generalization of the density functional theory [52-61] for relevant observables such as the electromagnetic current 
[62, 63] and the electromagnetic field. The current and electromagnetic field profiles, given by the equations of motion, 
help us to perform a partial resummation of the perturbation expansion. The study of the expectation values at scales 
£ Ri should give us some new insight into the quantum-classical transition because on the one hand, the description 
includes a large number of degrees of freedom needed to cope with classical objects and on the other hand, it is based 
on expectation values which are the relevant quantities at both sides of the transition. Notice that this is not the case 
when the usual, scattering amplitude based formalism of Quantum Field Theory is used because it is restricted to 
pure states and the decoherence can not even be formulated properly. Finally, we see macroscopic physics for ^„ < £. 

Three aspects of microscopic classical theory will be mentioned in this Section. One is the decoherence and another, 
induced by it, the transfer of the time arrow by the environment. The third issue is the treatment of the polarization 
effects of the Dirac-see. 



A. Decoherence of non-relativistic particles 

Let us consider, as a simple but non-trivial problem, a system of A'' non-relativistic charges interacting with the 
photon field. The dynamics of charges is characterized by the action Sc[x] and the coupling to the electromagnetic 
field is realized via the current 

N 

jI^U = (122) 

n=l 

where x^^^t) denotes the trajectory of the n-th charge. We consider the functional integral 

Z = J £)[i]£)[^]ef ^•-Do"'-^+»s[^^]-ieiWa-A ^^23) 

where D[x] denotes the integral over the trajectories of all charg ) stands for the particle trajectory, 

S\x] = /S'[a;+] — 5'*[a;~] and the current j[x] = {j[x~^],j[x~]) couples to the photon field. 

The photons can be eliminated in the usual way and the result is the action at a distance form of the electromagnetic 

interactions 



'[£]e'^'«[*l, (124) 
with 

S,s[x] = S[x] + W'P''°'[-ej[x+],ej[x-]] (125) 

where the last term on the right hand side is the influence functional [15]. The form (28) for the generating functional 
W^OBcb] gives 

soBc[x] = s[x] - ^f[x] ■ or ■ r[x] - yr w • . nx] - ^mi^Dr'' + iQDo)r[x] (ue) 

for open boundary condition where the parameterization (14) of the current coupled by the matrix a gives j'^ = 

and j'^ = — k)/2 + + k)/2. The relations (30), (40) allow us to rewrite the exponent of the integrand for 

open boundary condition in terms of the effective actions as 

5oBC [x] = S[x]+ rgj,"* [A, A^''] + irgj,"* ™ [^™] - ef [x] ■ {A + a™) - e'f [x] ■ A'^^ , (127) 
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where the fields A, A'^'^^ and A™ are the solutions of the classical equations of motion (33), (42). The classical field 
theory based on the effective actions rP B°* [A, A""^^] and 3°* [A™] allows us to obtain the Wilsonian effective action 
for the charges in terms of the expectation value of the electromagnetic field, A, and other non-physical expectation 
values, A'^'^^ and A™, induced by the trajectories x. 

The imaginary part of the influence functional (126) is due to the continuous spectrum. It results from the 
construction of the reduced density matrix, the measurement of the unit operator, and represents a suppression 

g = g40[»=+l-j[»="l)-3^o-0[x+]-j[a!-]) (^28) 

of the off-diagonal fluctuations in the density matrix p(a;+, x~). The imaginary part of the two current point function 
receive contributions form mass shell intermediate states, real photons, in this case. Non-accelerating charges do not 
emit photons; it is easy to check that s = 1 in this case. 

Let us now consider a charge moving along a the trajectory x{t) = rcosojt, the corresponding current being 
jt,x — ~ 3^(i))(l) 'i(f)). We shall calculate the imaginary part of the influence functional for jf-.,. = jt,x±R/2, 
r = \r\ <C \R\ = R'm Fourier space where 

(^eB{q\io,qr),r'-^[Biq°-u,,iv,q-r)-B{q°+uJ,u,q-r)]^ (129) 

with 

B(^^w,a) = ^e-'«''*+*«=°^'^* = (5(9°)-Fm(5(9°-w)-Fia(5(g°-Fw)-FO(a2), (130) 
and the suppression factor can be written as 

The diffraction integral in the exponent shows clearly the origin of the decoherence, the interference between the 
photons emitted by 7+ and . The power scries on the right hand side of Eq. (130) corresponds to the multipole 
expansion and we keep the dipole field only. Straightforward calculation gives in the limit tf — U ^ 00 

when R and r are parallel. For R - r = one finds 

= g4e=(t/-ti)r=[-^sinwr-^coswr-4] ^ ^^33^ 

The suppression factor interpolates smoothly between sy = s_l = 1 for i? = and 

S|| =si=e-t^'(*/-**)'-'-' (134) 

for ii pa oo. 



B. Radiation time arrow 



The Lagrangian of QED is formally invariant with respect to the inversion of the direction of time, nevertheless our 
daily experience confirms that the currents generate retarded electromagnetic fields and the time reversal invariance 
is lost [64]. This symmetry breaking must come either from the environment, represented here as external sources 
which couple to the current and the electromagnetic field or from the way the dynamical problem is posed ie. the 
boundary conditions in time. 

The CTP formalism is well suited for the investigation of this problem because of the duplication of the degrees of 
freedom. In the decomposition (7) of the field variable cj) picks up the diagonal quantum fluctuations in the density 
matrix which appear in the expectation values of the canonical coordinates and (j)'^'^'^ collects off-diagonal fluctuations, 
displayed by the expectation value of the canonical momentum variables. 

But there is another way of looking at this separation: from the point of view of the time arrow. Let us return to 
the remark made in Section II that the time arrow is different for the coordinates 0^ appearing in the time ordered 
and anti-time ordered part in Eq. (3). Complications arise from the fact that matrix elements of time dependent 
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operators include both advanced and retarded effects. For instance, the Green function {0\T[(p{x)(j){y)]\0) displays 
both advanced and retarded parts. The discussion of the linear response formulae in Section II makes it clear that any 
reflecting final boundary conditition induces advanced fields. Can one separate the retarded and advanced effects? 
The answer is trivial and affirmative for non-interacting particles and the decomposition (7) gives and (j)'^^ as 
the combinations of the canonical coordinates which receive retarded or advanced effects of the external sources, 
respectively, for open boundary condition, according to Eqs. (35) and (36). The retarded nature of the physical field, 
(j), can be understood by noting that the in-field, observed at the initial time, is vanishing due to the choice I*,) = |0). 
The separation of the advanced and retarded effects is trivial because it can be achieved by inspecting the initial (in) 
or final (out) field in the absence of vertices. 

Interactions mix advanced and retarded effects because the non-linear pieces of the effective action in the fields 
are time reversal invariant. Consider for instance the equations of motion arising from the effective action r[J, ^] of 
Eq. (117). The time reversal invariance is broken by the environmental variables a and j which appear only in the 
terms linear in the fields, the part of the effective action which depends on the boundary conditions in time. The 
relation between the expectation values J and A, given by the boundary condition independent quadratic part of the 
effective action, involves the time reversal invariant near field propagator. The local expectation values appear to 
be insensitive to the symmetry breaking effects of the boundary condition. Our conjecture is that closed quantum 
systems maintain time reversal invariance for any boundary condition in time. 

Such a claim makes the coupling to the classical environment responsible for the choice of the time arrow. But the 
environment characterized by the external sources a and j obey quantum laws, as well, suggesting that it is the classical 
limit for the environment which triggers the breakdown of the time reversal invariance. In fact, the decoherence which 
is supposed to be the hallmark of the quantum-classical crossover suppresses advanced fields without influencing 
retarded combinations, according to Eq. (126). Notice the double role the imaginary part of the influence functional 
plays in Eqs. (126)-(127). It sets not only the direction of time of the decoherence but determines the proportion of 
retarded and advanced parts in the charge propagator, too. 

The equations of motion of the effective actions T[J, J^'^^ , A, A'^'^^] support the view that the retarded and advanced 
propagators appear in the equations of motion due to decoherence. On the flrst sight it is rather confusing that the 
quadratic part of the effective action in the fields involve the retarded and advanced Green functions when open 
boundary condition is used, in contrast to r[J, How can different equations of motion occur at all for the same 
expectation values? The only difference between the equations of motion for J and A, arising from the effective 
actions r[J, J^'^ , A, A'^'^"] and r[J, ^] is that different quantities are kept fixed during the variations. In fact, let us 
construct a variation of A by varying the external sources a and j for fixed a and j in such a manner that J remains 
unchanged. Such a variation gives rise to the equation of motion for A, Eq. (120), of the effective action r[J, 
In the case of the effective action r[J, J'^^, A, A'^'^] all four sources are varied by keeping all field but A fixed. The 
lesson is that the time arrow appeared in the equations of motion relating measured expectation values because the 
variations respected the condition J'""^^' = A^"^^ = 0, the vanishing of the off-diagonal fluctuations. This is exactly 
what decoherence is for. 

Can the decoherence solve the arrow of time problem? Not, it makes one step only in reducing it to the environment. 
When the photons are eliminated in Section VI A then their boundary condition in time fixes the sign of S-Dq and 
consequently the arrow of time for the charges. What is remarkable is that this transfer of symmetry breaking becomes 
enhanced as we enter the macroscopic domain by decreasing the cutoff. This suggests that a better understanding of 
the arrow of time problem can be obtained by the application of the renormalization group method to this problem. 



C. Vacuum polarization 

In classical electrodynamics the polarizability of a medium is taken into account in the well known, simple manner 
supposing the higher multipole moments are negligible. It is interesting to see what kind of approximation corresponds 
to this procedure when the polarizability of the Dirac-see is considered. 

The effective action for the electromagnetic field, 

T[A] = ^A ■ {D'^''^'-^ - e^GS*^^"-^ + e'^G^''^'-^ ■ W^j^^^'^"'"'' ■ G^''^'-^) -A-eA- r^"^ + 0{A^) + O(e^), (135) 

obtained by eliminating the current J by its equation of motion and reintroducing the elementary charge e in the 
minimal coupling in Eq. (117). We ignore the gauge fixing term because current conservation decouples the gauge 
sector and write r[A] = T[F], F^j^i, = d^A^ — d^A^j^ due to gauge invariance. A conserved current is defined as 

lS_m -ifs'^^-a^m^ (136) 

' \ / |F=0 
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The minimal coupling is constructed by the introduction of the tensor H/j^i, = d^ji, — d^j^ for the current and writing 
the 0{A) part of the effective action as 

-eA.j = ^F^,-D-'-H^'-' (137) 

where the conservation of the currents was used on the right hand side. The field strength tensor which includes the 
vacuum polarization is defined as 

To find the linearized equation of motion for the photon field it is sufficient to have the accuracy 0{J) and O(A^) for 
the effective action, 

r[A] = -^F-G-eA-j + OiA^). (139) 

To find a more detailed expression for the field strength tensor G return to the effective action T[A] and assume 
that its 0{A^) part is translational invariant. This latter property allows us to write 

1^ . T . r(2) .t.A = -If^^. ^[r(2), F]f - If^, -L^- ^[T^^\Fr + If,^ -L^- ^[r(^), F]r (140) 

where 

([r, F]±)r = / [r[^%F^ ± i^^r, (i4i) 

•'y 

the integration over the space-time coordinate y being shown explicitly. We write F^^^ = —Dq^'^^~^+Y., cf. Eq. (135), 
S being the photon self energy, and find 

G = F + ^p, - \l . [E, U-^F]_ + \{L- [S, U-^F].r (142) 

where the superscript tr denotes the exchange of the two space-time indices in the last term and renders G an 
antisymmetric tensor. 

ft is instructive to obtain the field strength tensor in three-dimensional notations when Af^ = {(pj—A), E = 
—da A — V^, B = V X A and the parameterization 

E,. = (^^° ^) (143) 



yield 



^ ~ "^2 l^(So-n)£;-s X B S(g)£;-£;os- [n,$[jB]]+ 

,1^-2/ d{2do^-E + d-v)-do{dou + dw) . 

2 \-d{2doT. ■E + d-v) + do{dou + dw) dod'^u^ - dod^u^ + d^d^w^i - d^d^w^^) ^ ' 



with 



u = -(n + So)^; - S X B 
V = -(n-FSo)-E-l-S X B 

w = S(g)£;-i-£;0S- [n,$[B]]_ (145) 

f^[B]jk = ijkiBe, being the spacelike part of the bare field strength tensor F. 

The motion of the charges in the Dirac sec mixes the electric and the magnetic fields in the vacuum polarization 
if translational invariance is assumed only. The rather involved form of the polarization shown in Eq. (144) results 
from the Schwinger-Dyson resummation of the one-loop photon self energy parameterized as (143). It is interesting to 
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note that similar resummation has already been used in classical electrodynamics in arriving at the Clausius-Mossotti 
formula [65]. 

The polarization simplifies considerably when rotational invariance, Ujk = ^5jk, S = imposed. Expression (144) 
results in 



D = E + n-^ (^Eo - • [1 + n-\d^ - a d)]^ e 



and 



H = B + U- 



(So - s„)B + -s„,n-iao9 X E 



(146) 



(147) 



respectively, where D and H are the electric and the magnetic fields contained by the field strength tensor G, and 
Enr = n + Sq. 

Further simplification is achieved by imposing Lorcntz invariance. Enr = 0, because the electric and magnetic 
sectors no longer mix. The dielectric tensor e and magnetic permeability tensor ji defined by the relations D = e ■ E 
and B = ^ - H, respectively, are 



-1 = i + n-iSo, 



(148) 



reproducing well known one-loop result [66]. Notice that the operator D^^ on the right hand side is the result of 
gauge invariance, the relation between the gauge field and the field strength tensor. 

The n-th multipole moments of the vacuum polarization arc constructed by projecting the £ = n angular momentum 
components out of the current n-point Green fmiction. The dipole contribution corresponds to the retaining of the 
O(a^) part of the generating functional iy[a, j]. The order of the retained dependence on j gives the accuracy in 
terms of the powers of e, according to Eq. (C3). The effective action formalism where the Legendre transformation 
performs the inversion of the quadratic part, the resummation of the Schwingcr- Dyson equation, yields the relation 
between the applied and the polarized fields. The total field strength (144) is based on the dipole contribution to the 
polarization taken into account in the same order in e as the photon self energy (143). Higher multipole moments 
introduce higher order terms in the fields in the effective action and lead to non-linear polarization effects. Note that 
the n-th order multipole contributions do not appear before the order n of the expansion in e. 

We have traced the effects of vacuum polarization in the real part of the influence functional for the photon field 
which give the quantum corrections to the classical Maxwell equations. The imaginary part of this functional eg. SJe 
in a relativistically invariant theory plays important role, as well, by introducing decoherence and setting the time 
arrow for the electromagnetic fleld. 



VII. QUANTUM RENORMALIZATION GROUP 

The problem of the quantum-classical transition can be studied by adjusting the resolution of the space-time averages 
which is usually realized within the framework of the renormalization group method. This method was developed 
in its full generalities first when the Kadanoff- Wilson blocking was applied to the description of critical phenomena 
in classical Statistical Physics [67]. The application of blocking in Quantum Field Theory followed immediately and 
agreement with the already known multiplicative renormalization group scheme was established. 

However, there is an essential difference between the successive elimination of degrees of freedom in a classical 
and a quantum system which was ignored. In classical statistical physics the form of the partition function can be 
maintained during the blocking but the vacuum-to-vacuum transition amplitude of a quantum system is changed in 
a fundamental manner when particle modes are eliminated. In fact, the transition amplitude between two states 
corresponds to a process occuring for pure states while in Quantum Mechanics the elimination of a degree of freedom 
generates mixed states to be described by a density matrix. This goes beyond the usual formalism based on "in-out" 
type vacuum-to-vacuum transition amplitudes. 

Let us follow the blocking procedure on the path integral (8) where some high energy particle modes, p > A, are 
integrated over. As long as the the initial and final states are pure the density matrices factorize in the variables 
0^ and the successive integration over particle modes of the fields (j)^ preserves this factorization. Notice that such 
a blocking corresponds to the elimination of the particle modes which possess given initial and final states coded by 
the density matrices. Let us suppose, for example, that the occupation numbers for the high energy modes to be 
eliminated are fixed in the initial and final states. Then the resulting path integral displays a dynamics for the low 
energy particles which depends on these occupation numbers. 
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This is not what is meant by chmination of a degrees of freedom in Quantum Mechanics. The measurement of 
the identity operator, the calculation of the trace of the density matrix in the factor space of the degrees of freedom 
to be eliminated, destroys all information about the initial and final states. The information about the final state 
is destroyed by the trace operation, the dependence on the initial state disappears due to the unitarity of the time 
evolution. 

It is this point, the realization of such a true elimination of degrees of freedom without retaining information about 
the initial and final states, where the CTP formalism becomes essential. We erase the final state information for the 
high energy modes to be eliminated by installing the identity operator as density matrix in that sector of the Fock 
space. The identity operator, the functional Dirac-delta for the given modes of — (f>~ couples the variables along 
the two time axis and this coupling is handed over for the low energy modes. In other words, once we perform a true 
trace operation on the high energy sector of the Fock space the remaining low energy states become mixed. This 
mixing requires the use of the CTP formalism for the rest of the blocking procedure. 

The difference between the blocking in the vacuum-to-vacuum amplitude and in the CTP path integral is that in 
the former case the trace operation is performed in the zero particle number sector of the high energy modes only. 
This represents a small error when the the particle modes to be eliminated are at very high energy compared with 
the observational scale and perturbation expansion applies. But in non-asymtotically free models where the high 
energy dynamics is non-perturbative or in effective theories with not too high cutoff A the weight of the particles at 
the cutoff scale is non-negligible and the difference between these two blocking becomes more important. It is the 
quantum renormalization group procedure, based on the CTP formalism, what should be followed in these cases. 

There is no difficulty in extending the fimctional renormalization group scheme [68, 69] to the CTP formalism 
[39, 40] but the doubling of the degrees of freedom renders the resulting renormalized trajectory formal. A more 
natural formalism, based on the measurable expectation values and the degrees of freedom controlling decoherence, is 
outlined in Appendix E. But instead of embarking a detailed study of this evolution equation we shall be satisfied by 
a few qualitative remarks about the emergence of classical physics. The renormalized trajectory of realistic models 
passes at several fixed points and may experience different scaling laws separated by crossover regimes [70] . One of 
the crossovers belongs to the quantum-classical transition. The difference between the trajectories generated by the 
blocking in the traditional, scattering amplitude based formalism and the CTP path integral should be less important 
at the UV side of the quantum-classical crossover. The decoherence which requires the density matrix is supposed 
to be the hallmark of the quantum-classical crossover. Therefore one expects that this transition and possible others 
further down towards the infrared direction are missed completely or reproduced with substantial errors by the "in- 
out" formalism which can not handle mixed states. Note that we need a bosonizcd form of the theory where the 
excitations are handled by bosonic fields in order to embark the problem of the quantum-classical crossover because 
the expectation values of fermionic operators are trivially vanishing. 

The decoherence suppresses the fluctuations of the field (f)'^'^^ . But it is just the first moment, the expectation 
value of this kind of fields which is suppressed by construction in deriving the variational equations of motion of 
the effective action r[J, J^"^^ , A, A'^'^^]. This suggests that these equations which already assume decoherence change 
less at the quantum-classical crossover than the variational equations of the effective action r[J, A] which do involve 
non- vanishing advanced fields. 

It is worth mentioning that the O(A^) part of the action r[J, A], given by Eq. (117), has already been suggested 
many years ago [48-50]. By analogy of the steps leading from Eq. (100) to Eq. (105) we find for vanishing external 
sources the well known action at a distance 

Kr[J] = -- J • G"^''"'-! • J = r'"<^'='^[J] + ^J- -Do''''' • J, (149) 

with time reversal invariant current-current interaction term. The earlier proposals were made on the basis of classical 
electrodynamics where the charges follow world lines and the self-interaction was excluded in order to avoid UV 
divergences. In Eq. (149) self-interaction is included. Such theories have no radiation field, accelerating charges 
induce near field only. The retarded and advanced fields appear with the same weight and therefore the action- 
reaction ballance is obviously satisfied. The well established retarded radiation field can be recovered by taking into 
account the reaction of all charges to the electromagnetic field in a completely absorbing Universe [71] . As mentioned 
above, the boimdary conditions in time influence the linear part of the action only. As a result, the 0{A'^) part is the 
same for open ended time evolution and for fixed boundary condition when the photons end up at the vacuum state 
in a completely absorbing Universe. In other words, the classical argument of Ref. [71] applies to open boundary 
condition imposed on the quantum level, as well. The time arrow and the Abraham-Lorentz force [46, 72-74] are 
generated by the presence of other than the accelerating charge in the system. But this argument is valid in the 
microscopic domain only. As soon as the cutoff is lowered and our action is given in the classical domain, decoherence 
is implied automatically and the time arrow is generated by the retarded Lienard-Wicchert potential of the action 
r[J, J^'^", A, A^^^^j and the Abraham-Lorentz force is generated by the accelerating charge itself. 
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The renormalization group approach reveals a similarity between spontaneous symmetry breaking and the arrow of 
time problem. Spontaneous symmetry breaking appears in the renormalized trajectories as a crossover to an infrared 
scaling regime which contains relevant operators of lower symmetry. The boundary condition in time imposed on 
photons is passed over the charges when photons starts to behave classically at the quantum-classical crossover. At 
the final count, the arrow of time problem appears as a spontaneous breakdown of the time inversion symmetry. 
As such, it indicates that an extreme sensitivity develops in the system among the infrared modes to choose a time 
arrow from "outside". The identity of thermodynamical arrows formed in different domains and the importance of 
the cosmological master arrow can be explained by such a mechanism. 

We close this Section by pointing out an unexpected difficulty arising at the quantum-classical crossover. When 
reaching the crossover by decreasing the cutoff the advanced coordinates become suppressed. The semiclassical limit 
corresponds to the dominance of the path integral by the configurations (f)'^'^ = 0. But there is still another variable, 
(j>, to integrate over. Notice that there is a cancellation between the actions 5'[(/)"'"] and S'*[(/i~] in the exponent on the 
right hand side of Eq. (8) when 0*^"^^ = causing the integrand of the path integral to become flat. This appears in 
the renormalization group scheme as a strong coupling regime for the configurations (j). Thus, decoherence induces 
strong coupling behavior and our approximations obtained within the framework of the loop-expansion break down. 
Even the degrees of freedom may change. 

This scenario is in strong contrast to the semiclassical limit discussed usually in the framework of the path integral 
formalism. The limit ^ of the single time axis non-relativistic path integral 



is treated by the stationary phase approximation arguing that the integral is dominated by the paths close to the 
classical trajectory. This description applies to a closed quantum system which is not realistic in the macroscopic 
limit. When a density matrix is introduced to describe the separation of the system from its environment, the time 
axis is doubled and the CTP path integral 



is used where Sinfi[4>~^ ,4'~] denotes the influence functional of the environment. This expression is not dominated 
anymore by the classical trajectory. Instead there is an approximate cancellation between the two time axis when 
decoherence sets in. 5*1,111 ~ M^((t>^ — (6^)^/2 with large and the physical field, cf), is integrated over with approx- 
imately flat integrand. The classical limit is supposed to be recovered on the level of the equations of motion for the 
expectation values of observables only. The approximation to the transition amplitude (150) in the semiclassical limit 
is based on the rigidity of the dominant trajectory, the excitations being suppressed by 1/h. In contrast of this simple 
picture the semiclassical limit is supposed to be reached when the energy level density becomes high. This view takes 
into account the environment and is supported by the strong coupling, soft saddle point scenario of the functional 
integral (151) of the CTP formalism. 



The dynamics of the expectation value of local operators in QED, the electric current and the electromagnetic field 
in particular, was studied in this work. The expectation values were given in the framework of the CTP formalism 
covering either fixed initial and final states or open, unconstrained time evolution. 

It was shown that the equations of motion for the expectation values of local operators can be obtained as variational 
equations of a suitably defined effective action mixing the diagonal and off-diagonal quantum fluctuations in field 
diagonal basis. The effective action plays a double role. It is not only a classical action for a classical field theory 
for the expectation values of local operators but it appears also as a Wilsonian effective action for degrees of freedom 
coupled to the system by our observables. 

The effective actions for free photons and electrons are calculated first. While the former is a trivial exercise the latter 
proves to be a highly involved problem due to the multi-particle aspects of the Dirac-see. The dynamics of localized 
charged states depends on the states which contribute to the current expectation value in a manner reminiscent of 
small and large polarons of Solid State Physics. On the one hand, negative energy single-particle states contributing 
to the polarization generated by weak external field overlap strongly. As a result of the dynamics of two-particle 
states generated by the current from the vacuum, Pauli-blocking generates a long-range correlation which produces 
separation-independent forces preventing the separation of localized charges by weak external electromagnetic field 
in the non-interacting Dirac-see. On the other hand, the charges arising from localized states within the mass gap 




(150) 




(151) 
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decouple from the filled up negative energy states and from each other. For such charges the standard relativistic 
action for localized particles was recovered. The actual value of the mass parameter depends on the details of the 
bound state and it gives rise a renormalization condition on the classical field theory level. 

The effective action for the interacting electron-photon system was constructed up to quadratic parts in the fields 
on the two-loop order. The elimination of either the current or the electromagnetic field produces the one-loop action 
at a distance theory. The quadratic part of the one-loop effective action for the current and the electromagnetic field 
calculated in the leading order of the gradient expansion agrees with the action of classical electrodynamics. Quantum 
corrections to the action come from three directions. Vacuum polarizations due to electron-photon interactions 
generate form-factors and new, higher order contributions in the fields which can be treated pcrturbativcly having 
as small parameter. Other e-independent quantum corrections appear due to the indistinguishability of electrons in 
the non-interacting Dirac-see. The small parameter to organize an expansion in the exchange effects is the absolute 
magnitude of the ratio of the electron Compton wavelength square and the invariant length square of the separation 
in space-time. Finally, the third class of contributions corresponds to the boundary conditions in time which are build 
in at a microscopic level. 

The Dirac-see is the natural polarizable medium of QED and its effects can be taken into account in the same manner 
as in classical electrodynamics, by the introduction of the deplacement fields and dispersive electric permeability, when 
the action of the classical field theory is truncated in the quadratic approximation. 

Though the effective actions mentioned in this work govern expectation values of a single measurement only they 
give some qualitative informations as far as multiple non-demolishing measurements are concerned. The correlation 
function of two successive non-demolishing measurements of the current is the inverse of the O(J^) kernel of the 
effective action. In general, the multiple current measurement results are weakly correlated by vacuum-polarization 
effects when they are carried out in space-time locations separated by an invariant lengths whose absolut magnitude 
exceeds the Compton wavelength. 

The CTP formalism is particularly well suited to trace the way the quantum boundary conditions in time generate 
the time arrow in classical physics. A conjecture was put forward, namely that closed quantum systems maintain time 
reversal invariance for any boundary conditions in time. When the system becomes open by coupling it to external 
sources the response is given in terms of retarded Green functions and a time arrow is generated for unconstrained, 
open ended time evolution. It was found that advanced effects mix in at some order only if the system is not allowed 
to follow an open ended time evolution and is projected on some states at the final time. Advanced effects represent 
an influence of the future on the present and are possible in Quantum Mechanics due to the superposition principle 
which allows us to specify both initial and final states. The symmetrical treatment of two time axis with opposite 
sense of time in the CTP highlights the fact that such effects are as natural as the influence of the initial condition 
on the course of the motion. Similar phenomenon exists in classical physics, as well. By specifying the initial and 
final coordinates of a classical system the change of the final coordinate modifies the whole time evolution. But 
this is a triviality because the state of the system became uniquely defined at the final point only, the equations of 
motion being of second order. Notice that the influence of the boundary conditions in time does not decrease when 
the distance in time between an observation and the boundary conditions is increased, in other words there is no 
clusterisation in time due the unitarity of the time evolution. 

The CTP formalism duplicates each degree of freedom, (t> ^ (j)^, for the treatment of the density matrix p{(p'^,(j)~) 
and the two copies are subject to opposite time arrow. The sum -|- ^~ (difference t/)"*" — </>") remains invariant 
(changes sign) under the inversion of the time and follows retarded (advanced) dynamics. The observed averages 
are given by a certain linear combination of and 0^ and classical field theories constructed for the measured 
fields only can not identify the sense of time. As a result, the dynamics of such theories is time reversal invariant 
by mixing retarded and advanced effects with equal weight. The retarded and advanced effects can be separated 
by retaining both t/)"*" and in the classical dynamics. The resulting dynamics reflects decohcrcncc on the average 
because (\E'|0^ — 4>~\'^) = 0. The classical field theory constructed for both and (p~ predicts indeed retarded 
dynamics for the measured combination. Decohcrcncc, the suppression of <j)^ — <j)~ , generates macroscopical time 
arrow by suppressing degrees of freedom subject to advanced dynamics. 

The effective action formalism offers classical dynamics both for the microscopic and macroscopic regimes because 
the space-time resolution of the expectation values is limited by the UV cutoff A. For high enough A we sec the 
microscopic quantum structure of local operators and the effective action is an extension of the density functional 
theory. As the cutoff is lowered one recovers a hydrodynamics description, classical field theory for the current, 
and macroscopic physics. The issue of quantum-classical crossover can be addressed in this generalization of the 
renormalization group scheme for the density matrix. The decoherence is specially clear in the CTP formalism where 
it appears as a simple one-loop effect. It is conjectured that the quantum-classical crossover which is characterized 
by decohcrcncc is actually a strong coupling phenomenon for the physical fields. 

The qualitative similarity between suppression of advanced combination of the degrees of freedom and decoherence 
suggests that the equations of motions of classical field theories which trace both degrees of freedom of the CTP 
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formalism changes less at the qiiantum-classical crossover than those of the classical theories containing the measurable 
fields only. Thus the equations of motion for the measurable field interpolate between the time symmetrical form of 
the classical theory for the physical fields only and the retarded equations of the classical theory for all fields as A is 
decreased and cuts through the quantum-classical crossover. This reasoning shows the peculiarity of the Abraham- 
Lorentz force of radiation damping. This force originates from the consistent motion of all charges, interacting by the 
near-field Lienard-Wiechert potentials when the space-time resolution of the expectation values is microscopic. But 
once the space-time resolution is chosen to be macroscopic the Abraham-Lorentz force arises from the accelerating 
"blocked" charge alone. Such a view of radiation damping applies to dissipative forces in general and explains the need 
of auxiliary variables when such effects are incorporated into canonical dynamics. Another remarkable fact is that 
the quadratic part of the action for the physical fields is independent of the boundary conditions in time. Consider 
two Universes corresponding to the same initial conditions but one is allowed to have an open, unconstrained time 
evolution while the state of the other is projected on the vacuum at a final time. The latter represents a completely 
absorbing Universe and its linearized equations of motions are identical with that of the Universe with open ended 
time evolutiuon. 

The application of Quantum Field Theory to establish the dynamics of the expectation values offers a comprehensive 
framework for measurement theory. The basic difficulties of the measurement theory is to identify well established and 
understood quantum effects in the interaction of a small and a large system. The difficulty comes from the presence 
of a dimensionless number, n/N, the ratio of the number of degrees of freedom in the two systems taking extremely 
small values and producing unexpected effects in a manner analogous to Statistical Physics. The n/A/'-dependence 
of the expectation values can be traced by the study of the cutofF-depcndence because the small system can be 
imagined as residing in the elementary volume element of the underlying field theory model. The cutoff-dependence 
can systematically be obtained by eliminating the degrees of freedom, by the application of the renormalization 
group strategy based on a suitable chosen blocking procedure. The quantum-classical crossover shows some formal 
similarity with spontaneous symmetry breaking, both representing new, less symmetrical scaling laws in the infrared 
scaling regime. The radiational time arrow is passed between subsystems by this mechanism. In the case of QED 
the time reversal invariance for the charges is broken by the boundary conditions in time imposed for photons at the 
quantum-classical crossover. This mechanism is realized by the imaginary part of the photon propagator which is 
non-vanishing for continuous spectrum, ie. the infinite time limit only. Another, conjectured feature of the quantum- 
classical crossover, the strongly coupled dynamics, might be part of the problem which prevents us from having a 
better insight into the quantum-classical transition. 



This work could not have been carried out without the stimulating discussions with Janos Hajdu whose critical 
reading of the manuscript was helpful, too. Furthermore I thank Kernel Sailer for the careful check of some formulae 
and Andor Frenkel for pointing out the relevance of refs. [13, 14] . 



In this Appendix the generating functional is constructed for a hermitian scalar subject to boundary conditions 
Pi = |0)(0| and p/ = 1. We start by collecting some formulae for the propagators of any hermitian, spinless, local 
operator (f). The retarded, advanced near and far field Green functions are introduced in the next step, followed by 
a brief recall of the spectral representation. Finally, the propagator of a free, neutral scalar particle is discussed and 
the massless results are given in closed form. 
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APPENDIX A: TWO-POINT FUNCTION FOR HERMITIAN SCALAR FIELD 



1. Propagator as a block- matrix 



The connected Green functions are obtained by means of the generating functional 




(Al) 



and the propagator is defined as a 2 x 2 block-matrix 



A"" = - 

xy 



(A2) 
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When both legs of the propagator are on the positive-sense time axis one finds the usual causal propagator 

iA++ = J2{n\T[<l>M0)mn) 

n 

= {0\T[<p,cPy\0). (A3) 
The propagator along the negative-sense time axis is 

iA-- = ^(n||0)(0|f[^,</),]|n) 

n 

= {Q\T[cf>ycf>,]\Oy 

= (iA++)*. (A4) 

Finally, the mixed propagator is given by 

^Ai'y" =E<0|'^^l^)H</'x|0). (A5) 

n 

The choice p/ = 1 allows us to decrease tf until it reaches max(x°, and placing both field operators into the same 
expectation value resulting 

iA+- = {0\<f>y4>^\0). (A6) 
Note that the identity Tlcp^^iy] + T[(j)x(f)y] = (px't'y + 't'y'Px leads to the relation 

A++ + A— = A+- + A-+. (A7) 

2. Retarded and advanced propagators 

Introducing the notation A = A++ for the causal propagator 

iA,,,. = @{t - t')(0|.^,</.,nO) + e(t' - i)(0|.^x'<^x|0), (A8) 

its real an imaginary parts are 

JRA,,,, = -ie(t- f')«(O|[0x,<^x']|O) - ie(t' - t)i(0| </.,]|0) = -ie(f- i')»(0|[</'x,</'.']|0) (A9) 
with e{t) = sign(t) and 

9A,,,, = -ie(i-O(O|{</'x,'^x'}|O>-^0(i'-O(O|{'/'c.',<^x}|O) = -i(O|{0„<^,4|O), (AlO) 
respectively. The retarded and advanced propagators are defined as usual, 

iAt, = ±e{±{t-t')){o\[ct>.,Mm- (All) 

It is easy to find the actual expressions in terms of A by noting 

iA"- = e(f-O(o|0x0x'|o> + e(f-t')(o|0x'0x|o) 

= lA,,, - e{t' - t){o\<^x>ci)x\o) - eit - i')(o|</'x'</'x|o) 

= iA^a!>-iA+-, 

= iAl^,+iA-+, (A12) 

where the relation iA = (iA++)* and Eq. (A7) were used in obtaining the last equation. The near and far field 
propagators. A" = KA and A-'' = 25iA satisfy the equation 

A'' = A" + iA^. (A13) 
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Similar reasoning yields the equation 

A° = A" - ^A^ (A14) 

for the advanced propagator. The complete propagator can be parameterized by means of three real functions, for 
instance 

^ - ViA/ + iQA -A" + iQA ) ^^^^^ 

3. Spectral representation 

A particularly useful parameterization of the propagator is provided by the spectral representation. By following 
the standard procedure introduce the eigenvectors of the energy-momenta, P"|n) = Pn\n, to induce the space-time 
dependence (px = e'^^(^oe~'^^ and write 

i{0\(Px(l^y\0) = ^i(0|(/)^|n)(n|</.^.|0) 

n 
n 

= i j e-'f(^-^)e(/)p(p2), (A16) 
Jp 

where 

e(p°)p(p^) = ^^(27r)^<5(p„-p)|(0|</.o|n)r, (A17) 

n 

is the spectral function. One can recast the last line of Eqs. (A16) in the form reminiscent of free propagators, 

Jp Jo 

= i [ I dM^e-''^''(''''(*-*''+'''('=-"'V(M'), (A18) 

Jk Jo 



where u)k{ii^) = -\//?+^ and 

dk 



(A19) 
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The propagator A is obtained by means of the Fourier representation of the Heavyside-function, 

Axx, = - / d^i^p{^l^) / — + — . (A20) 

and performing the change of variable lu —uj in the second integral 

rco p /-i(ujk(fj.^)+uj){t-t')+ik{x-x') gi(wfc(;i^)-w)(t-t')+ifc(a;-x') \ 
A^^' = - / / —. + 



(J + ie —u) + ie 



where 



/>OQ 

/ dMW)Aoxx'(M') (A21) 
Jo 

The steps above repeated for the block-propagator yield 

POO 

A= / d/xV(M')Ao(M'). (A23) 
Jo 
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4. Free neutral scalar particle 

The path integral representation of non-interacting generating functional is 

= j D[^]ei^-^^'-*+''^-^. (A24) 
where Aq = Ao(m^). The path integral can be carried out with the result 

Wo\j\ = ~3-^o-l (A25) 
The detailed form of the propagators can be obtained by means of the Fourier integral 

0^= /[e-''^'=*+''=^afc + e^'^'=*-^'="'4] (A26) 

Jk 

where ujk = lOkijn^)- The non-vanishing canonical commutation relation are [ak, a\,] = {2'n)^2ujk5k,k> ■ The repetition 
of the steps (A20)-(A21) leads to the causal propagator (A22). 

5. Real space expressions for m = 
First we calculate the expectation value 



k 



1 

— dkk d cos 616-''=*+''='' '^""^ (A27) 

Jo J-1 



87r2 

The contributions to the /c-integration are negligible at the upper limit when working in a cutoff theory which can be 
incorporated by calculating the integral fov t ^ t — ie yielding 

(010.^.10) = (A28) 

and 

Do = Ao(0) being the free massless propagator. The retarded and advanced propagators are 



D^o = ^0(±(^° - ^°'))^((^ - ^0')- (A30) 



APPENDIX B: TWO-POINT VERTEX FUNCTION FOR NON-INTERACTING DIRAC-SEE 

In this Appendix we present the calculation of the kernel of the quadratic effective action for the current with /? = 
for the non-interacting Dirac-see. 

The Fourier representation of the quadratic form (71) 

rSi,,,,.„,^-^ie(«.-y + ete-wiX. -■'<-■--■> l^j^^^ (bi) 

gives after straightforward integration over the frequency and the polar angles 

e(i)r(2)-i(t,r) = + r) - zt/i(t - r) + h{t + r) - h{t - r)] (B2) 
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with t = ti — t2, r = \xi — X2\ and 



fn{t) 



-I- 



ikt 



(B3) 



The relation 9"/„(t) = {—i)^5{t) yields fi{t) = Ci — iQ{t) and /2(t) = —t@{t)+tC2 + C3, Ck being integration constants 
to be determined. The vertex function 



r(^h\t,r) = i^^^[2c2 - e{t - r) - e{t + r)] 



(B4) 



obtained by replacing the solution for and /2(i) into Eq. (B2) should be symmetric in t therefore c = 1/2 and 
Eq. (72) follows. 

APPENDIX C: LOOP EXPANSION FOR W[a,j] 

In this Appendix some details of the calculation of the generating functional for connected Green functions of the 
current and the photon field in QED are presented. 



1. Setting up the loop expansion 
We start by Eq. (87) for the generating functional written after the photon field is integrated out, 



The identity 



yields 



where 



6f U/ ^/ ; 



e 2 



m=0 



fj-Do-j 

^ ml 



a' = eaDo ■ j 



(CI) 



(C2) 



(C3) 



(C4) 



and 5' /5a' acts only on the first factor of W^^. For simplicity the counterterms are suppressed. The j-dependence 
comes from a! therefore we have 



m=0 



1. 



(C5) 



2. Two-loop order quadratic generating functional 



We shall calculate W[a, j] on the two-loop order up to quadratic terms in the external sources. To this end we need 

W = W'ia + a'] + A . ie^aDoa ■ -^rrW^ + a'] + l(^- le^aDoa ■ ^) W% + a'] (C6) 
oa oa 2 \da da J 

and its square in Eq. (C5) up to two-loop terms. We start with the expansion in a', 

W^'ia + a'] « W^'ia] + W^'^'^ + i<f ^ [a]a'A + [aja'Ad', + l^'/j [aja'^a'^d', (C7) 



4! 
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by retaining the 0(a*) pieces only. Two functional derivations give 

^W'^'ia + a'] ^ W:'^'^ [a] + [aK + ^Cli'^ [^KK + ICli'i 

^W^^[a + a'] « <P[a] + <;f[a]a'e + iw^„f/i[a]a;a'., (C8) 
where the coefficients are 





[a] . 






^el(2) 

ab 


[a] . 






^el(3) 

abc 


[a] . 








[a] . 




, ■'■11/61(6) - - 



(C9) 



in the given order. The photon exchange appears in the combinations 



^el(4)| 
abed I 



W^^^'\aU = W:l%]{azDoaU. (CIO) 



The expressions 



i/ • '^Doa ■ + a'] « [a] + • a' + ^a' ■ ^^^^ [a] • a' 



(^±.ai2Doa.^yW'^'[a + a'] « (<f [a] + <11^ (aiDoa)„e^. (Cll) 



give 



+e^W'^^^^ [a] + e^W'^^^^ ■ a' + ^a' ■ W'^^*^ [a] ■ a' (C12) 
in the given truncation. Its square 

^ (w%] + Wi'^[a]a', + ^W^\a]a'X^ + e\d + a') ■ W^''\a] ■ aiD^a ■ W^''\a] ■ {a + a') (C13) 

is obtained by taking into account unitarity, W'^^fO] = 0. We have now everything needed to re-exponentiate the 
necessary terms in Eq. (C5), 



l + iW - ^W'^ « expi 



W[a] + W^'^ [a] • a' + ia • W^^^ [a] ■ a' + e^W'^^^^ [a] + e^W'^^^^ ■ a' + y a' • W'^^*^ [a] ■ a' 



yielding 



W[a,j] = ~j-Do-j-^{a + a') • (5o • (a + a') + Y^abcd"'cO'd{^iDo&)ab 



(C14) 



+^'^^abcd"'a^d{^^^0^)'>c + Y^abcd i^^^O^)('daa"'b " — + a')aGoabi^Doa)beGoefia + d')f.{C15) 
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The substitution of al given by Eq. (C4) results in 

W[a,j] = ■ [-Do - e'^Doa ■ (Go - e^W^f^) • aDo - e^Do^ ■ Go • aDoa ■ Gq • aDo] ■ j 



(C16) 



-2a(-Go + e^W^^^ - e^Go ■ aDoa ■ Go)d + ej[-Doa ■ (Gq - e^W^^^) - e^Doa ■ Go • aDoa ■ Go]a 



leading to expressions (88)-(92) for the generating functional. The parameterization (14) and (67) give 



mV[a,a,j,j] = -{a,j,a,j) ■ i^^ ^tr g 




with 



-K'oBC = — 



Qret 
^0 



eG'"'^Dl''^ 



(C17) 



(C18) 



for OBC shown in a detailed form in Eq. (95). For FBC the sources belonging to different time axis decouple giving 
Eq. (C17) with 

_ _ / G"*^"*^ e(G'^«^"£)S<^^'' + 3G3£»o)' 

~ Ve(£>g^^^G"^'^^ + QDo^G) Do 

For the sake of completeness we record the imaginary part, 

e(9GL'g^^ + G^'=*9£»o)\ . fa"' 



(C19) 



and 



QWfbc 



1 / SjG 

^WoBc[a,j] = --{a,j) . (^^(^.etcjG + ^DoG^) 

( 



[a,a,j,j] = --ia,a,j,j)- 



SGI ^ ^ 
K 2 



e(£)^^'*^SG + QDoG"""^') ( « 



QD 



e(3GDS°''' + G"'="9i:'o 

K 2 



(C20) 



2 
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a 



i 

J \J. 
(C21) 



APPENDIX D: LEGENDRE TRANSFORMATION 

The calculation of the quadratic effective action is sketched in this Appendix. There is a useful relation between 

the second derivatives of a function and its Lcgendre transform. Let us denote all sources and fields in the problem 
by j and 0. Then the effective action is defined by r[0] = W[j] — j ■ (p with </> = 6W\j]/Sj. One can easily prove the 
identity 



6^W\j] 6^r[ 



SjSj 54>5(j) 



= -1 



(Dl) 



which renders the calculation of the Lcgendre transform of a quadratic expression to an inversion problem. 

We start with a few useful relations for the inversion of block-matrices to be used later and rules for obtaining the 
retarded and advanced components from a product of propagators. This is followed by the calculation of the effective 
action for four fields with open and fixed boundary conditions. 



1. Inversion of block-propagators 



We shall meet two kinds of inversion during the Legendre transformation of different functionals. One of them is 
the inversion of a 2 x 2 block propagator 



K 



G eSx 
eS2 D 



(D2) 
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The inverse 



-eS2 -D / V-eD™*-iS'2G-i 



(D3) 



52 -JJ J \-eU"" ^t>2( 

where 

= G-e^SiD-'^S2 (D4) 

and 

£)mt ^^j-e^SaG-i^i (D5) 
is obtained by solving the system of linear equations 

i) = ^ {i} 

by elimination. 

Anther inversion problem one encounters involves the block-matrix 

«<^^r ' . (D7, 

The same strategy as in the previous case ie. the solution of the corresponding set of linear equations gives 

2. Light-cone as a null-space 

The real part of massless particle propagators are non-vanishing on the; liglit-c;onc' only. This property simpli- 
fies enormously the expressions of the perturbation expansion frequently involving a chain of products of massless 
propagators 11 = Dia ■ D2(t ■ ■ ■ Z)„ and its retarded and advanced parts. 

We start with the demonstration of the useful relation 

£,far . ^near-l ^ q ^g-) 

valid for massless propagators. The straightforward proof of Eq. (D9) is by inserting the propagator 

ret p p — ik{x — x') 

with F{0) = 1 resulting 

(£,far . ^near-l)^^, ^ _ Pe-»'=(^-^') (DH) 



which is indeed vanishing for e — > 0. 

To understand better the rcgularization effects of e we note that on the one hand Dq = (Z)"^*** -|- D^'^)/2 satisfies 
the inhomogeneous equation of motion 

^ • □ • Do = 1- (D12) 



F{~U) 



Thus, G^'^^^ oo and Gnear-i ^ Q ^Yie light cone for e ^ 0. On the other hand, D^'^" is non-vanishing on the light 
cone only. Therefore, Eq. (D9) follows. A more explicit argument is based on the fact that Z)^*'' is a homogeneous 
Green function satisfying the equation 



^ . □ . £)far = (D13) 



F{-a) 
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which is just Eq. (D9) due to Eq. (D12). 

The lesson is the identity Dn>=ar-i . jjret ^ ^^near-i . ^^adv^ •g_ ^.jjg equation D'"'*' = D'^" holds when the propagators 
appear in a product beside of a factor £)"ear-i rj,j^-g j-ggyj^ yields 

^near— 1 ^ ^^et ^^near— 1 ^ j^adv ^ 

Let us finally consider the product 11 where each propagator is of the form 
in the basis (j+ — + It is easy to verify that 11 is of this form, as well, and 

n,A...n - ( Dfv.2?adv...^adv 

(o ^Di ■ Df^ ■ ■ ■ Df^ + Df* ■ ■ ■ ■ Df^ ■■■+ Dl^* ■ Df* ■ ■ ■ Qd) ^^^^^ 

giving 

{Dia---DnT = Df^---D^^. (D17) 

The property (D14) can be used to show 

^ (Df * • • • D'^* + Df^ ■ ■ ■ Df^) = D^^''' ■ ■ ■ D^'"'' (D18) 

giving 

{Dia ■ ■ ■ DnT^''" = -Or"' • • • Dl""'- (D19) 

The fact that the product Dia ■ ■ ■ D„ preserves the form (D15) assures that W^'^''^ which appears in the self energy 
is of this form, too. 

Eqs. (D17) and (D19) can be used to show that once the pair creation is excluded (9G = 0) the generating 
functional WFBc[S,i] is obtained from WoBc[a)i] by the simple rule found for non-interacting systems, namely the 
replacement of all retarded or advanced propagator by their near field version. 



3. Effective actions for four fields 

To obtain the effective action for the fields J, J'^'^, A and A'^'^^ we need the inverse of the block-matrices of Eqs. 
(95) and (96). We start with OBC when the two-loop expression electron self energy in Eq. (91) and the first equation 
in (D17) give 

grret ^ Qrct ^ ^2Qrct . ^rct . Qret _ g2^el(4)ret_ ^20) 

In a similar manner the photon self energy of Eq. (92) gives 
the inverse being 

= Dj;''*-! - e2(Q*'* - W'^'^'^^"''). (D22) 

These expressions substituted in Eqs. (D4) and (D5) give D'"* = D'^^^ and G'"' = Gj^''* - e'^W'^^^^"'\ The inverse 
(D3) gives the expression (111) when choosing 5*1 = (5''°*" • £'0''* = ^2- Finally, the effective action (110) follows from 
Eqs. (Dl) and (D8) and its equations of motion give 

Gga = -{l + e'^W^^^^''-G^-^)-J'^'' + eG^-A'^^, 
Df^j = eDf^-J'"^^-A'^^, 

or ■ (« + f «) = -(1 + e^W^'^"^-^^* • Gf • (j - ^r"-) + ■ {a ^A^'^') , 

DT ■ (j + ^j) = ■ {j - f -Z'^'^) -A+ (D23) 
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According to the remark made at the end of Appendix D2 the effective action rFBc[>^j'^ j^, ^ ] is obtained 
from FfbcI'^) J^^^, A, A'^^] by the replacement of all retarded or advanced propagators by their near field version as 
long as pair creation is neglected. 



EfFective actions for two fields 



The effective action which involves physical fields results from the Legendre transformation in the variables a and 
j only. The generating functional ^W[a,a,j,j] based on the two-point function (95) or (96) is 



_ ^near 



eG 



near _ jjncar 



Wa 



-D 



.ret 



■J 



Do J \Jy 

where the only boundary condition dependent pieces are 

WoBCa = -G'"' ■ {a + eZ?f • j) WoBCj = 
or 

WpBCa = -G'^"^'' • (a + eD^^'^'^ ■ j) WfeCj =-Do-j- eZ^g^^-^ • G"<^^^ • a. 
The application of the inverse (D3) produces the effective action 



g^ret . Qret 



rFnc[J,A] = -^iJ,A)-K-'- 



l(J,A).if-i 



Wa 

Wj 



with 



K- 



1 fGT'"'-^ + G^""'-^ ■ w;:'" ■ G, 



(4)n /Vnear-1 
D ■ '^0 



— e 



— e 

T-jnear— 1 
-^0 



(D24) 

(D25) 
(D26) 

(D27) 
(D28) 



APPENDIX E: FUNCTIONAL RENORMALIZATION GROUP IN THE CTP FORMALISM 



In the original version of the rcnormalization group method one follows the evolution of coupling constants one- 
by-one in the framework of the perturbation expansion [67] and exact results might be achieved by resumming the 
perturbation expansion. The functional realization of the renormalization group [68, 69] is based on one-loop evolution 
equation and the exact results should arise by letting infinitely many term mixing in the action. In this Appendix 
we outline briefly the evolution equation for the effective action in the latter formalism which can easily be converted 
into a numerical, non-perturbative algorithm to solve models. 

Let us introduce a cutoff k in the theory for the field variable (f> in the quadratic part of the action, S'[(/)] — » 
S[(j)] + i(j) ■ Kk ■ (p/2. The properties of the real operator K arc: (i) It should suppress all fluctuations for k ^ oo, ie. 
the eigenvalues of K^o should be oo. (ii) The physical theory should be recovered for A: = which is guaranteed by 
the condition Kq = 0. (iii) For finite, non-vanishing k K/. should suppress the modes with momentum below k ie. 
the eigenvalues of the translation invariant Kh should be large for momenta p < k. The generating functional for the 
connected Green functions of the field is 



I 



D[4 



,iS[<l>]-^,l,-K-4,+ij-^ 



(El) 



(j) — ((/)+, (j> ) and S[(p] = S[4i'^] — 5* [0 ]. The generating functional is equipped with UV regulator which is omitted 
for simplicity. The quadratic term in the exponent controls the fluctuations and 



K 



K 
K 



The evolution of the generating functional 



W = -Ti 
2 



Sj+dj+ Sj Sj- 



s^w .5\vsw_ .m^sw_ ^ 

Sj+ Sj+ Sj- Sj 



(E2) 



(E3) 
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is obtained by bringing the derivation with respect to k, denoted by a dot, into the functional integral. The parame- 
terization (14) of the external sources gives 



W = Tr 



fS^W 



6^W 1 



.-2 X2 



We need the evolution of the effective action 

r[<^,<^n=r[, 

where we separate the real and imaginary parts 



SjSj 



)\-k 



xadvi 



im j,adv im] 



^W[j,j] 



ladv 



3 ■< 



ladv im 



(E4) 



(E5) 



(E6) 



and introduce the variables 



xadv 



ladv im 



SW\j,j] 

6j ■ 



(E7) 



The sources are real and therefore, two field variables are independent only. We shall use (j), (j)^'^ as independent 
variables and introduce the functional T^™[(f>,(j>^^] = Qr[(j>^™,(j)'^'^^ ™] for the imaginary part. The derivation of the 
evolution equation for the effective actions is straightforward except that wc have to keep track of the evolution of the 
relations among the independent and the dependent field variables. The expressions which follow will appear simpler 
by the help of the notations $ = (0, 0*^"^), (^™ = (?!>™, 0*"^"^ ™) adopted in the rest of this Appendix. 
The equation (Dl) allows us to rewrite the evolution equation (E4) in terms of the effective action, 



Tr 



+i 



Kl 



4 \S4>S$), 



(0+0+ + 



jadv im-Aadv im 



(E8) 



■K 



were the subscript indicates block-matrix elements of the second functional derivative. Furthermore the disconnected 

contributions to the two-point fimctions arc given in terms oi cf>^ = cf) + ?'0™ — (k =F 1)(0^'*'^ + icf)'^'^^ ™)/2. 

The evolution equation (E8) is not yet closed because it contains the dependent field variables, too. To relate 0'™ 
and locally in the field configuration space we introduce the derivative matrix 



Sab[4'] = 



50; 



im 



The elimination of the variables 0™ is achieved by the help of the derivation of the equations of motion 

. (53?r 6T 



(E9) 



(ElO) 



with respect to 



■ s. 



The derivation of this equation with respect to A yields the evolution equation for S, 



Sab 



(XL)- 



'd6 • 



cd 



de 



'^eb T^K^ Jeb 



(Ell) 



(E12) 
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After the solution of the problem posed by the dependent variable we return to the evolution equation (E8) and 
we separate the tree-level contribution by the parameterization 



r = f + iSjf + -(0+ ■ k ■ + 4)- ■ k ■ 4)-) 



(E13) 



of the effective action. The real and imaginary parts of V satisfy the evolution equations 

-1 




+ 



S + B 



1 + k2 / ^2sjjp 



4 
-1 



A 



+ ■ 



■K 



■S + B 



^adv im^adv 



(E14) 



■K 



respectively, where 



Bab = 2K„.h — 2 




-K. 



cm 



(E15) 



The effective action of the physical system is obtained by integrating the system of equations (E12), (E14) from 
fc w 00 where a perturbative initial condition is imposed down to A; = 0. The renormalized trajectory displays the 
scale dependence of the effective action. In the case of QED one should start with jjroperly regulated theory eg. the 
effective action should be calculated as in Eq. (82) after the replacement D'^^ Dq^ + Tk\ in the action for the 
initial value of A but keeping the original, A-independent counterterms. The projection T into the transverse photon 
states in the suppression term renders the turning on of the interactions of photons with the fermion loops according 
to their momentum in a gauge invariant manner. 
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